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abstract 


A  theoretical  method  is  derived  whereby  the  effect  of  the  compressible 
chree-dimensional  relief  on  the  torque  required  for  a  helicopter  rotor 
jay  be  calculated.  The  complementary  wing  approach  is  used  to  represent 
the  relief  effect.  The  effect  of  the  wing  which  complements  a  finite  wing 
to  make  it  an  infinite  wing  may  be  found  as  a  reduction  in  the  free-stream 
velocity.  This  approach  was  applied  to  fixed  wings  by  Anderson1  and  is 
extended  here  for  application  to  helicopter  rotors.  By  using  a  simplified 
blade  element  analysis,  the  magnitude  of  the  change  in  torque  due  to  tip 
relief  has  been  calculated  and  compares  well  with  test  results.  The  method 
to  use  the  tip  relief  calculation  in  more  complex  blade  element  computer 
programs  is  described. 


PrBCBdiRS  Pa6e 


ui 


FOREWORD 


This  investigation  was  sponsored  by  the  Eustis  Directorate,  U.S.  Army  Air 
Mobility  Research  and  Development  Laboratory,  under  Contract  DAAJ02-71-C- 
0006,  Task  1F162204AA4201 ,  during  the  period  from  October  1970  to  October 
1971. 


The  basic  concept  underlying  this  investigation  was  suggested  by  Mr.  James 
P.  Trant,  Jr.,  Eustis  Directorate,  USAAMRDL,  whose  continued  assistance 
during  the  program  is  gratefully  acknowledged.  Appreciation  is  expressed 
to  Dr.  Gabriel  D.  Boehler  for  the  guidance  provided  during  the  investigation. 
Also  acknowledged  are  the  detailed  notes  and  program  decks  made  available  by 
Mr.  Robert  Sopher.  Further,  a  limited  distribution  paper**1  was  made 
available  by  the  Douglas  Aircraft  Company,  McDonnell  Douglas  Corporation. 

The  technical  representative  of  the  Contracting  Officer  for  this  contract 
was  Mr.  Paul  Mirick. 


Preceding  page  blank 


v 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  . . .  ill 

FOREWORD  .  v 

LIST  OF  ILLUSTRATIONS  .  ix 

LIST  OF  TABLES  .  X 

LIST  OF  SYMBOLS  .  xi 

INTRODUCTION  .  1 

BACKGROUND  .  3 

Rotary-Wing  Performance  Calculations  . 3 

Compressibility  Effects  on  Rotors  . . .  4 

The  Three-Dimensional  Relief  Effect  on  Rotors  . .  6 

THE  TIP  RELIEF  EFFECT  .  8 

Physical  and  Mathematical  Models  for  the  Tip  Relief  Effect  .  8 

The  Correction  for  a  Semi-Infinite  Wing  . . . . .  11 

The  Correction  for  a  Finite  Wing  . . . . .  18 

Airfoil  Section  Characteristics  . . 20 

Characteristics  of  the  Correction  Factor  .  23 

Correction  Factor  for  a  Rotating  Wing  . . . .  31 

Correction  for  the  Change  in  Dynamic  Pressure . . .  32 

BASIC  EQUATIONS  OF  A  ROTOR  IN  COMPRESSIBLE  FLOW  .  37 

Fluid  Dynamic  Equations  for  a  Rotor .  37 

Direct  Solution  . 43 

Prandtl-Glauert  Type  Transformation  for  a  Rotating  Blade  .  46 

Extensions  of  the  Compressible  Flow  Transformations  .  49 

APPLICATION  OF  THE  TIP  RELIEF  EFFECT  IN  ROTOR  PERFORMANCE 

CALCULATIONS  .  50 

The  Physical  Concept  of  Tip  Relief  Applied  to  a  Rotor  Blade 

Element  Analysis  . 50 

Tip  Relief  Calculation  in  Blade  Element  Analysis  .  53 

A  Simplified  Blade  Element  Analysis  .  54 

Vi.  i 


Preceding  page  blank 


RESULTS 


Page 

58 


CONCLUSIONS  .  63 

LITERATURE  CITED  .  64 

DISTRIBUTION  .  70 


viii 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Coordinate  System  for  a  Finite  Wing  .  9 

2  Coordinate  System  for  a  Se-ai-Infinite  Wing  . .  11 

3  Geometrical  Relationship  . 22 

4  Coefficients  of  the  Taylor  Series  for  a  Semi-Infinite 

Wing  .  25 

5  Coefficients  of  the  Taylor  Series  for  Finite  Wings  of 
Three  Aspect  Ratios  With  a  10%  Thick  Parabolic  Arc 

Airfoil  Section  . 26 

6  Correction  Factor  for  a  Semi-infinite  Wing  With  a 

Parabolic  Arc  Airfoil  Section;  Thickness  ■  10%  .  27 

7  Correction  Factor  for  Finite  Wings  of  Three  Aspect 

Ratios  With  a  10%  Thick  Parabolic  Arc  Airfoil  Section  .  28 

8  Effect  of  Aspect  Ratio  on  the  Correction  Factor 

Calculated  by  the  First  Term  of  the  Taylor  Series  .  29 

9  Distortion  of  a  Stream  Tube  Over  a  Finite  and  an 

Infinite  Wing . . . 33 

10  Theoretical  Correction  of  the  Drag  Coefficient  vs. 

Mach  Number  for  an  NACA  0012  Airfoil  From  Two- 

Dimensional  Data  to  an  Aspect  Ratio  of  3  .  36 

11  Coordinate  System  .  42 

12  Prandtl-Glauert  Type  Transformation  of  a  Rotor  .  48 

13  Transformation  of  a  Rotor  Blade  Element  .  52 

14  Change  in  Torque  Coefficient  Due  to  Tip  Relief  for  a 

Hovering  Rotor  (a/b  ■  0.0159) . . . .  60 

15  Change  in  Torque  Coefficient  Due  to  Tip  Relief  for  a 

Hovering  Rotor  . 62 

ix 


Table 


LIST  OF  TABLES 


Page 

I  Comparison  of  Integral  Accuracy  With  Number  of  Terms  . .  14 

II  Airfoil  Constant,  In/t  . . 23 

III  Correction  Factor  for  a  Semi-Infinite  Wing  . .  30 

IV  Change  in  Torque  Due  to  Tip  Relief  .  59 


x 


list  of  symbols 


Defined  in  equation  (2),  area 

Speed  of  sound 

Defined  in  equation  (112) 

Number  of  blades 
Semichord 
Drag  coefficient 
Pressure  coefficient 
Drag,  distance 

z-coordinate  of  wing  or  blade  surface 
Correction  factor 
Defined  in  equation  (4) 

Defined  in  equation  (23a) 

Defined  in  equa~:  (19) 

Defined  in  equativ  (25) 

Defined  in  equation  (26) 

Source  of  strength 

Logarithm  to  the  base  e 

Mach  number 

Pressure 

Dynamic  pressure 

Total  radius  of  a  blade 

Radiol  station 

Nonaimensiorial  radial  station 
Reynolds  number 


t 

-* 

V 

vxlVy*vz 

WVz 

x,y,z 

6 

Y 

C.n,5 

X 

y 

6a 


p 

a 

T 

4> 

SUPERSCRIPTS 

O 
(  )' 

SUBSCRIPTS 

(  >x 
<  )« 

<  >i 


Sami  thickness 
Frea-stream  velocity 
Velocity  vector 

Velocity  components  in  Cartesian  coordinates 
Velocity  components  in  cylindrical  coordinates 
Coordinate  system  (see  Figure  1) 

Prandtl-Glauert  factor 
Ratio  of  specific  heats 
Auxiliary  coordinate  system 
Aspect  ratio 
Advance  ratio 

Ratio  of  area  of  airfoil  profile  and  of  enclosed  rectangle 

Ratio  of  moment  of  Inertia  of  airfoil  profile  and  of 
enclosed  rectangle 

Density 

Blade  solidity 

Thickness  ratio 

Potential  function,  inflow  angle 

Averaged 

Derivative,  auxiliary  coordinate,  nordimensional 

Differentiation  with  respect  to  x 
Condition  in  the  undisturbed  stream 


xii 


Incompressible 


Effective 


Finite  wing 
Semi-infinite  wing 

Differentiation  with  respect  to  the  auxiliary  variable  £ 
Chord,  compressible 
Refers  to  local  value 


xiti 


INTRODUCTION 


The  present  methods  of  calculating  the  performance  of  helicopter  rotors 
in  the  compressible  regime  are  not  quite  satisfactory.  These  methods  are 
based  on  the  use  of  two-dimensional  compressible  airfoil  data  inserted  into 
strip  element  analyses  and  predict  power  losses  higher  than  are  actually 
found  from  flight  test  data. 

Though  only  one  of  many  complex  factors,  a  prime  reason  for  this  discrepancy 
is  the  tip  relief  effect,  i.e.#  the  reduction  in  the  drag  of  a  finite  body 
when  compared  to  a  similar  body  of  infinite  span,  due  to  the  reduction  in 
local  velocity  near  the  tip.  The  existence  of  the  tip  relief  effect  has 
been  well  documented  for  both  fixed  and  rotary  wings,  and  the  documentation 
indicates  that  the  tip  relief  effect  becomes  more  important  as  Mach  number 
increases  in  the  subsonic  range.  For  the  case  of  a  wing  in  uniform  flow, 
theoretical  work  exists  which  quantifies  the  tip  relief  effect.  Further, 
a  theory  exists  by  which  the  drag  coefficients  of  wings  of  differing  aspect 
ratios,  up  to  high  subsonic  speeds,  may  be  correlated  by  a  similarity  law 
developed  by  Anderson.1 

The  purpose  of  the  work  reported  herein  was  to  investigate  whether  Anderson's 
work  on  the  tip  relief  effect  for  fixed  wings  could  be  extended  to  rotating 
wings,  either  by  a  transformation  of  coordinates  from  uniform  to  rotating 
flow,  or  by  setting  up  the  problem  directly  in  rotating  coordinates.  The 
approach  was  suggested  by  Trant2  several  years  ago.  The  nonlifting  finite 
wing  is  represented  as  a  source  distribution,  and  the  actual  wing  is  found 
tc  be  mathematically  equivalent  to  the  superposition  of  a  wing  of  infinite 
aspect  ratio  and  of  a  "complementary  ’.ring"  that  would  extend  from  the  tin 
of  the  actual  wing  to  infinity;  the  problem  is  initially  sot  up  for 
incompressible  flow  and  then  transformed  into  a  compressible  flow  problem. 

Early  in  the  investigation,  it  was  found,  exactly  as  theorized  earlier  by 
Trant,  that  Anderson's  method  could  be  used  for  rotating  wines  by  using  a 
different  method  of  averaging  locally  induced  velocities  over  the  wing  due 
to  the  tip  effect.  While  Anderson  finds  a  unique  correction  factor  for 
the  free-stream  velocity,  one  can  find  a  spanvise  varying  correction  factor, 
which  can  be  used  in  the  blade  element  performance  analysis  of  rotors.  Only 
recently,  however,  has  a  full  understanding  of  the  problem  been  reached, 
after  (1)  the  mathematical  transformation  of  the  problem  from  fixed  to 
rotating  coordinates  was  accomplished  and  (2)  a  solution  of  the  problem 
directly  in  rotating  coordinates  was  attempted.  The  limits  of  validity 
of  the  solution  indicated  above  have  been  established,  and  it  has  been 
shown  that  the  Trant  approach  is  correct  for  helicopter  rotors ;  in  addition, 
there  is  strong  evidence  that  the  mathematical  formulation  of  the  Anderson 
problem  directly  in  rotating  coordinates  is  froitless.  limited  comparisons 
of  flight  test  data  with  a  simplified  blade  element  analysis  incorporating 
the  tip  relief  effect  derived  theoretically  in  the  reoort  show  very 
encouraging  correlation  trends. 
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This  report  presents  a  detailed  mathematical  formulation  of  the  tip  relief 
effect  on  helicopter  rotors.  The  report  first  discusses  the  background  to 
the  problem |  then  the  calculation  of  the  spanwise  varying  correction  factor 
to  quantitatively  determine  the  three-dimensional  relief  effect  is  given. 
Whenever  possible,  the  discussion  and  the  derivation  are  ip  terms  of  fixed 
coordinates,  as  it  is  simpler  to  grasp  physically  and  because  the 
relationships  between  fixed  and  rotating  coordinates  are  straight  forward 
as  is  shown  later  in  this  report.  Next,  the  basic  fluid  dynamic  equations 
of  a  rotor  are  discussed,  including  the  derivation  of  the  small  perturbation 
equations  for  subsonic  and  transonic  flows  to  show  the  relation  between 
compressible  and  incompressible  flows  over  rotors.  The  source-sink 
distribution  method  of  representing  a  rotating  blade  ic  discussed.3  Then, 
use  of  the  tip  relief  correction  factor  for  compressible  flow  over  a  rotor 
in  blade  element  analysis  is  derived.  The  final  chanters  make  a  comparison 
vith  experiment,  discuss  the  results,  and  give  conclusions  and  recommenda¬ 
tions. 


BACKGROUND 


ROTARY-VINO  PERFORMANCE  CALCULATIONS 


Performance  prediction  techniques  for  rotary  wings  are  highly  iterative. 
Calculations  for  a  single  point  on  a  performance  curve  may  take  a  considera¬ 
ble  time  on  a  high-speed  computer  using  one  of  the  more  sophisticated  free- 
vake  techniques.  This  is  clearly  unacceptable  for  preliminary  design 
purposes,  and  thus  methods  which  are  reasonably  simple  in  application,  but 
not  necessarily  exact  in  the  mathematical  sense,  are  needed. 

The  two  approaches  to  rotsry-ving  performance  calculations  are  the  actuator 
disc  and  the  blade  element  theories.  Since  actuator  disc  theory  assumes 
an  infinite  number  of  blades,  it  gives  no  information  as  to  how  the  blades 
should  be  designed  and  ignores  profile  drag.  Thus  it  is  necessary  to  use 
blade  element  theory,  which  assumes  that  each  element  can  be  considered 
aerodynamieally  as  an  independent  two-dimensional  airfoil  segment.  The 
local  forces  are  calculated  from  the  resultant  velocities  at  the  element, 
and  total  forces  are  found  by  integrating  the  contribution  of  each  element 
along  the  radius  and  the  azimuth.  From  these  forces,  the  motion  of  the 
blade  is  determined,  but  the  motion  in  turn  has  an  effect  on  the  local 
velocities  and  thus  the  local  angles  of  attack.  This  iteration  may  be 
avoided  only  by  making  such  further  assumptions  as  averaging  the  induced 
velocities  and  assuming  a  constant  drag  coefficient  as  is  done  in  chapter  8 
of  the  book  by  Gessow  and  Myers.4  With  the  use  of  high-speed  computers,  the 
blade  element  iteration  is  used.  Some  of  the  basic  assumptions  used  in 
the  analysis  and  some  of  the  suggested  improvements  include  [for  a  more 
thorough  review,  see  Gessow  and  Myers4(page  198-216),  Harris  et  al.,5 
and  Bellinger6]* 

-  The  use  of  two-dimensional  steady  airfoil  data  is  almost  universal. 
The  use  of  so-called  synthesized  data  has  been  suggested  and  will  be 
discussed  later  in  this  report.  There  are  also  attempts  to  include 
unsteady  stall  hysterieis  in  the  data,6  in  addition  to  theoretical 
work  in  the  area  of  dynamic  stall.7 

-  Tip  effects  are  accounted  for  as  a  loss  in  the  lifting  capability 
of  a  small  percentage  of  the  blade  near  the  tip.  Often,  it  is 
simply  taken  as  a  constant  B»0.97,  but  it  can  be  calculated  as  a 
function  of  thrust  coefficient  based  on  vortex  theories  developed 
by  Prandtl  and  Betz,  as  explained  by  Glauert9  (page  261-269). 

-  The  radial  component  of  the  velocity  at  an  element  is  neglected,  and 
the  analysis  is  baaed  only  on  the  velocity  normal  to  the  blade  span. 
A  means  of  removing  this  limitation  has  been  suggested  by  Harris,5 
based  on  ideas  of  swept-wing  aerodynamics. 

-  The  induced  velocity  is  in  many  cases  assumed  constant.  The  major 
problem  is  the  prediction  of  the  location  of  the  trailing  vortices. 
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In  fixed-wing  aerodynamics,  the  trailing  vortices  are  assumed  to 
trail  aft  in  the  direction  of  the  free-stream  velocity.  On  a  rotor, 
the  trailing  vortices  are  in  more  or  less  the  shane  of  a  helical 
sheet,  due  to  rotation  and  to  the  nonuniform  lift  distribution. 

The  helical  shape  for  the  trailing  vortices  was  first  used  by 
Goldstein.9  But  the  helical  shape  gives  rise  to  a  great  deal  of 
interference  which  distorts  the  wake,  as  described,  for  example, 
by  Landgrebe, 1 0  * 1 1 

-  Compressibility  effects  are  taken  into  account  by  using  airfoil  data 
from  two-dimensional  compressible  flow  tests  as  discussed  below. 


The  complexity  of  the  rotating  wing  problem  is  such  that  only  recently  has 
any  real  progress  been  made  in  the  basic  theory.  In  fact,  Sear912  indicated 
that  the  blade  element  approach  is  valid.  Only  in  December  1970  was  it 
shown13  that  Sears  is  not  necessarily  valid  for  a  finite  span  rotating  blade. 
Boundary  layer  calculations  (such  as  that  by  Hicks  and  Naah14and  Clark  and 
Arnold i 1 s) use  the  two-dimensional  pressure  distribution  suggested  by  blade 
element  analysis  as  the  external  flow.  Unfortunately,  very  little  of  the 
basic  theoretical  analysis  is  of  direct  use  in  the  design  and  performance 
prediction  of  rotor  systems. 

COMPRESSIBILITY  EFFECTS  ON  ROTORS 


Some  of  the  approaches  to  account  for  the  compressible  flow  over  heliconter 
rotora  will  be  reviewed.  In  the  early  1950's,  two  basic  theoretical 
approaches  to  compressibility  were  made.  The  first  was  solutions  to  the 
actuator  disc  problem  in  compressible  flow. ife » 1 7 » 1 8  This  was  combined  with 
blade  element  theory  by  Laitone19  for  use  with  helicopters,  The  other 
approach  by  Head,20  was  the  use  of  the  Praadtl-Glauert  rules  for  the  lift 
coefficient  variation  with  Mach  number  in  a  simplified  blade  element  analysis. 
By  including  a  linear  decrease  in  lift  above  the  lift  divergence  Mach  number 
to  extend  the  validity,  the  ratio  of  compressible-to-incompressible  thrust 
coefficients  and  the  ratio  of  induced  and  total  torque  coefficients  for 
compressible  and  incompressible  flows  are  found.  Neither  of  these  approaches 
recieved  much  use.  The  actuator  disc  theory  vas  not  used  because  it  does 
not  give  any  details  concerning  the  airfoil  design;  Head's  report  received 
little  circulation;  and  the  use  of  two-dimensional  compressible  airfoil 
data  in  blade  element  analysis  vas  begun  at  the  same  time. 

Airfoil  data  obtained  in  two-dimensional  compressible  flow  in  blade  element 
analysis  were  used  to  find  the  effects  of  comoressibilitv  and  to  make 
corrections  to  the  available  incompressible  rotor  charts.  Typicallv, 
these  were  numerically  integrated  by  calculators  and  made  a  large  number 
of  assumptions.2 1 »22  The  culmination  of  this  effort  is  the  charts  and 
tables  prepared  bv  Tanner,23  29  using  a  highly  sophisticated  computer 
program  with  airfoil  data  from  "two-dimensional  wind  tunnel  tests  of  a 
production  blade  specimen  using  the  NACA  0012  airfoil  section". 
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The  inconbietency  of  ueing  tvo-dimeneional  data  in  rotating  airfoil 
problems  first  came  to  light  during  propeller  tetts.25  The  drag  divergence 
tip  Mach  number  of  a  propeller  vas  found  to  be  about  0.06  above  the  tvo- 
dimensional  result,  and  performance  teats  of  full-scale  helicopter  rotors 
in  hover  indicated  the  same  trend.26  Therefore,  Carpenter27*28  proposed 
the  use  of  " synthesised"  airfoil  data,  which  are  obtained  by  "adjusting  the 
airfoil  data  in  such  a  manner  that  the  calculated  performance  agrees 
closely  with  the  measured  performance".  He  stated  further  on  page  2  of 
Reference  28  that  "although  the  synthesised  data  are  derived  from  rotor 
hovering  performance,  it  can  be  reasonably  assumed  (until  proven  otherwise) 
that  the  data  are  equally  applicable  to  forward  flight  calculations". 

Charts  similar  to  Tanner's  were  prepared  by  Kisielowski  et  al.,  using 
the  synthesised  data  of  Reference  28,  but  these  latter  charts  appear  to  be 
rarely  used.  It  would  be  of  interest  to  compare  the  predicted  performance 
using  each  of  the  two  sets  of  charts  with  some  actual  flight  test  data. 

The  problem  of  helicopters  at  high  forward  speeds  is  that  compressibility 
effects  and  stall  occur  at  the  seme  time  on  the  rotor.  At  the  advancing 
tip,  the  forward  speed  and  the  rotational  speeds  add,  resulting  in  speeds 
of  approximately  Mach  one.  But  on  the  retreating  side  the  two  speeds 
subtract,  resulting  in  a  region  of  reverse  flow  and  areas  where  the  flapping 
blades  are  at  or  above  the  stall  angle  of  attack.  The  dividing  line 
between  stall  losses  and  compressibility  losses  is  not  indicated  by  test 
data,26  and  both  areas  must  be  taken  into  account  in  rotor  analysis. 

To  reduce  the  data  and  decrease  the  amount  of  flight  test  required, 
personnel  at  the  Edwards  AFB  Flight  Test  Center  have  been  developing  an 
approach,  mainly  empirical,  to  account  for  compressibility  effects.30- 
Drawing  on  fixed-wing  experience,  the  power  required  is  divided  into  two 
portions;  (l)  the  incompressible  poveri  and  (2)  an  increment  in  power  due 
to  compressible  flow,  which  depends  on  the  difference  between  the  critical 
Mach  number  of  the  blade  section  and  the  Mach  number  of  the  element.  The 
correlation  is  established  empirically  using  UH-1F  and  CH-U7A  flight  test 
data.30*31  Later  reports32*33  include  the  effects  of  power  losses  due  to 
retreating  blade  stall. 

Another  approach34  to  the  aerodynamic  problem  of  a  rotor  blade  in  compressi¬ 
ble  flow  has  been  taken  by  Bell  Helicopter  Company  in  its  design  of  the 
Model  309  King  Cobra.  The  rotor  system  embodies  an  asymmetrical  airfoil 
section  and  a  double  swept  tip.  The  airfoil  section  designed  by  Prof. 
Wortman35  is  based  on  the  aerodynamic  requirements  for  two-dimensional 
airfoils  found  from  basic  blade  element  theory.  Swept  wings  have  long  been 
used  to  delay  the  effects  of  compressibility  on  fixed-wing  aircraft. 

Similar  sweeping  of  the  leading  edge  of  the  blade  tip  will  reduce 
compressibility  effects  in  forward  flight.36 

Thin  blades  near  the  tips  to  reduce  the  drag  divergence  Mach  number  have 
also  been  considered  by  Bell.36  But  in  uniform  flows,  a  thin  airfoil  has 
a  lower  maximum  lift  coefficient  and  stalls  at  a  lower  angle  of  attack, 
which  in  a  rotor  would  increase  the  power  on  the  retreating  side.  The 
Sikorsky  S-67  Blackhawk  also  incorporates  the  sweeping  and  a  thinner 
section  of  the  blade  tip  to  reduce  compressible  flow  effects. 
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Surprisingly ,  though,  full-scale  vind  tunnel  tests37  indicate  that  the 
lifting  capability  of  a  thin-tipped  bladed  rotor  is  not  compromised,  while 
it  requires  significantly  less  power  at  high  advancing  tip  Mach  numbers  than 
a  blade  with  constant  thickness, 

THE  THREE-DIMENSIONAL  RELIEF  EFFECT  ON  ROTORS 


This  report  is  concerned  with  one  aspect  of  the  rotor  compressibility 
problem.  As  noted  above,  two-dimensional  airfoil  data  overestimate  the 
compressibility  losses.  One  possible  factor  nay  be  the  reduction  in  drag 
coefficient  of  a  finite  wing  in  compressible  flow  compared  to  a  similar 
two-dimensional  airfoil.  As  will  be  described,  this  relief  effect  may  be 
determined  in  a  quantitative  manner  and  applied  to  the  performance 
calculation  of  rotary  wings. 

In  a  three-dimensional  flow,  the  stream  may  deviate  from  the  main  stream 
in  both  directions  perpendicular  to  the  body  chord  line.  At  the  same  time, 
in  two-dimensional  flow,  only  deviation  in  the  direction  perpendicular  to 
both  the  chord  and  the  span  is  possible.  Therefore,  lesser  disturbances 
are  produced  by  three-dimensional  bodies  than  by  two-dimensional  bodies. 

(See  the  discussion  of  Figure  10.)  This  reduction  in  disturbance  velocity 
is  also  a  reduction  in  Mach  number,  and  thus  a  reduction  in  critical  Mach 
number. 

That  this  disturbance  reduction  in  a  three-dimensional  flow  increases  with 
free-stream  Mach  number  was  found38  analytically  for  flow  past  ellipsoids 
of  various  aspect  ratios;  the  difference  between  the  maximum  velocity  of 
an  elliptic  cylinder  (an  ellipsoid  of  infinite  aspect  ratio)  and  a  finite 
aspect  ratio  ellipsoid  was  found  to  increase  with  Mach  number  (see,  for 
example,  Figure  13-2  of  Shapiro's  hook,39)  On  airplane  wings  of  normal 
aspect  ratio  at  incompressible  speeds,  the  wing  is  considered  to  be  two- 
dimensional.  The  effect  of  aspect  ratio  is  found  from  the  Prandtl  lifting 
line  theory  as  an  additional  angle  of  attack  and  an  induced  drag  due  to 
the  effect  of  the  trailing  vortices.  The  induced  drag  of  the  wing  is  found 
to  vary  inversely  as  the  aspect,  ratio, 

Goethert140  extended  the  Prandtl-Glavert  transformation  between  compressible 
and  incompressible  flows  to  finite  wings.  The  result  is  that,  instead  of 
the  simple  affine  transformation  of  a  two-dimensional  body  to  reduce  it  to 
an  equivalent  incompressible  flew,  for  a  three-dimensional  body  all  lateral 
dimensions  must  be  reduced  by  the  amount  8  ■  J  1  -  M*  ,  For  a  more 
detailed  discussion,  see  Shapiro,39  article  13-2,  and  Jones  and  Cohen. *♦  1 
pages  b8-52. 

But  at  higher  speeds,  not  only  the  lift  but  the  drag  of  the  ving  is  affected 
by  aspect  ratio,  and  the  Goethert  extension  does  not  account  for  the  three- 
dimensional  relief  effect.  Therefore,  a  means  to  quantitatively  predict 
the  relief  effect  is  needed.  Anderson1  developed  a  theory  for  the  relief 
effect  on  fixed  wings  in  the  nonlifting  case.  This  theory  is  extended  here 
and  is  applied  to  rotary  wings. 
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An.-er.3on' s  theory  for  the  tip  relief  used.  &  "complementary  ving"  approach, 
the  complementary  ving  being  a  fictitious  ving  extending  from  the  actual 
ving  tip  to  infinity  on  either  side.  His  basic  proposition  is  that  the 
velocity  at  the  location  of  a  finite  ving  can  be  considered  to  be  the  same 
as  for  the  infinite  ving,  reduced  by  the  amount  which  would  have  been 
contributed  by  portions  of  the  infinite  wing  that  were  removed  from  the  flow 
(the  complementary  wings).  These  values  can  be  calculated  by  linearised, 
subsonic  compressible  flow  theory,  using  a  source-sink  distribution  to 
represent  the  thickness  effect  of  a  symmetrical  wing.  In  the  case  of  the 
fixed  ving  considered  by  Anderson,  the  effect  of  the  complementary  ving  can 
be  averaged  over  the  planform  area  of  the  finite  ving,  resulting  in  a 
"correction  factor"  which  is  a  function  of  aspect  ratio  but  not  of  position. 
The  application  of  this  correction  factor  is  "that  the  finite  wing  will  be 
considered  to  experience  the  same  static  pressure  to  total  pressure  ratio  at 
corresponding  points  as  the  infinite  wing,  but  only  when  measured  at  a  free- 
stream  velocity  which  is  greater,  by  the  correction  (factor)  calculated, 
them  the  free-stream  velocity  over  the  infinite  ving".  In  other  words,  we 
can  make  use  of  two-dimensional  airfoil  information  in  the  calculation  of 
three-dimensional  flows  through  the  use  of  the  correction  in  free-stream 
velocity.  This  will  be  explained  in  more  detail  in  the  chapters  which 
follow. 
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THE  TIP  RELIEF  EFFECT 


A  quantitative  means  to  calculate  the  tip  relief  effect,  so  that  it  may  be 
applied  to  rotary  vings,  is  derived  in  this  section.  The  derivation  will 
follow  the  concept  and  the  treatment  proposed  by  Anderson.1  The  major 
difference  between  Anderson's  reports  and  this  report  is  that  Anderson, 
being  interested  in  fixed  wings,  intended  to  obtain  the  correction  as  a 
function  of  aspect  ratio  X  alone,  while  here  it  is  necessary  to  obtain 
results  as  function  of  aspect  ratio  and  distance  from  the  tip,  ao  that  it 
may  be  applied  to  rotors. 

The  result  which  will  be  obtained  here  and  that  found  by  Anderson  is  that 
two-dimensional  results,  be  they  from  theory,  wind  tunnel  or  flight  test 
data,  may  be  used,  through  the  transformation  derived  here,  in  the 
calculation  of  three-dimensional  flows. 

The  discussion  here  vill  be  concerned  with  a  wing  in  a  subsonic  uniform 
airflow  with  a  constant  airfoil  section.  The  flow  is  assumed  to  be  sub- 
critical,  with  no  shock  waves,  and  inviscid  without  separation.  The 
airfoil  is  assumed  to  be  thin;  thus  the  thickness  of  the  airfoil  may  be 
represented  by  the  linear  theory  distribution  of  sources  and  sinks,  namely, 
that  the  source-sink  distribution  is  proportional  to  the  local  streamwise 
slope  of  the  airfoil.  Although  this  is  invalid  at  the  leading  and  trailing 
edges,  because  of  the  averaging  process  which  will  be  discussed  below,  no 
correction  needs  to  be  made.  The  discussion  is  further  restricted  to  the 
nonlifting  case,  i.e«,  the  effect  of  thickness  only.  It  is  believed  that 
the  major  cause  of  tip  relief  is  thickness,  rather  than  angle  of  attack 
or  camber. 

The  derivation  will  be  made  for  a  wing  in  uniform  flov,  as  the  physical 
characteristics  can  be  most  easily  explained  and  the  equations  are  more 
familiar.  At  the  end  of  this  section  it  will  be  shown  that  the  result  is 
the  same  for  the  rotating  case  as  may  be  expected  from  the  theory  developed 
by  Sears.12 

The  Goethert  extension  of  the  Prandtl-Glauert  rule  is  used  in  the  analysis 
for  application  to  fixed  wings.  The  factor  82  ■  1  -  M2  will  be  used  to 
transform  between  the  compressible  flov  and  the  equivalent  incompressible 
flow,  rather  than  any  of  the  higher  order  developments.  Extra  complexities 
of  the  higher  order  theories  cure  not  Justified  for  this  application.  Further 
extension  to  rotary  vings  la  discussed  in  the  next  chapter. 

PHYSICAL  AND  MATHEMATICAL  MODELS  FOR  THE  TIP  RELIEF  EFFECT 

Consider  a  symmetrical  airfoil  wing  in  inviscid,  incompressible  flov  at 
aero  angle  of  attack.  Due  to  the  principle  of  superposition,  the  velocity 
at  any  point  is  the  sum  of  the  free-stream  velocity  and  the  velocity  induced 
by  the  wing.  The  induced  velocity  of  a  finite  span  wing  may  be  represented 
by  a  source-sink  distribution  over  the  planform  area  of  the  wing.  But 
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^igure  1.  Coordinate  Cystem  for  a 
Finite  Ving. 


A  = 


because  of  superposition,  the 
finite  span  wing  may  also  be 
renresented  by  a  similar  infinite 
3pan  vin g  les3  "complementary 
winds'1.  The  complementary  wing 
is  a  similar  wing  extending  from 
the  tips  of  the  finite  wing  to 
infinity  on  both  sides.  The 
concept  of  the  complementary  wing 
is  illustrated  in  figure  1  for 
a  wing  of  chord  2c  and  suan  of 
2X.  9e  may  write  this  discussion 
in  equation  form  as 


0 

I  ■  0  x  t  /  U„A  dii  (l) 

*  -2X 

for  the  potential  of  a  finite 
vine  of  aspect  ratio  X,  with  flow 
along  the  x-axis  and  origin  at 
the  right  wing  tip,  where  A 
represents  the  9ource-sink 
distribution  integral 


+c 


Fc(0  dC 


J  (x-V)2  +  (y-n)2  +  z2' 


-c 


(2) 


where  the  subscript  t  indicates  differentiation.  This  coordinate  system 
was  chosen  so  that  it  may  later  correspond  to  a  helicopter  rotor. 


The  potential  of  the  finite  wing  using  a  pair  of  complementary 


d  3  + 


UooA  d'i 


wings  is 
dn 


(3) 


The  values  of  d  obtained  from  either  equation  (1)  or  (3)  have  to  be  iuon- 
tical. 

The  first  integral  in  equation  (3)  (with  limits  of  plus  and  minus  infinity) 
is  simply  the  potential  of  an  infinite  or  two-dimensional  airfoil,  and 
thus  can  be  represented  by  the  symbol  $ 

A  =oo 


The  last  two  integrals  give  values  which  vary  with  position  throughout  the 
flow  field.  If  an  everage  value  of  these  two  terms  can  be  obtained,  the 
problem  becomes  greatly  simplified.  Let 

.-2X 

C(X,x,y,z)  -  J  A  dr,  +J  A  jn  (1) 

-»  o 
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Then  equation  (3)  becomes 


<t>  =  U.x  +  -  lV?(X,x#y,*)  (5) 

The  velocity  in  the  streamwise  direction  may  be  pimply  written  from  equation 
(5)  as 


u  "  IT  =  U-  +  uX-»  "u«Gx  (6a) 


or 

u  *  Ujl  -  Gx)  +  uXa<B  (6b) 


which  may  also  be  written  as 


du 

U 


(6c) 


If  an  averaged  value  of  G  can  be  found  (denoted  by  7?  ),  then  the  velocity 
on  the  finite  wing  u  can  £e  compared  to  the  velocity  of  an  infinite  wing, 
but  at  the  new  free-stream  velocity  of  U^d-S  ).  Therefore,  the  finite 
wing  will  experience  the  same  static  pressureto  total  pressure  ratio  at 
corresponding  points  as  the  infinite  wing  at  a  free-stream  velocity,  which 
is  greater  than  the  free-stream  velocity  of  the  infinite  wing  by  the 
correction  fl  .  Therefore,  the  finite  and  the  infinite  wing  will  experience 
the  same  force  per  unit  span  when 


iui-ax) 


finite  U 
wing 


infinite 

wing 


(7) 


holds  for  the  free-stroam  velocity, 

The  correction  5  (X)  was  found  by  Anderson1  by  averaging  over  the  planform 
area  of  finite  wfng 

1  1  r  +c  r 0 

8  (X)  ■ - /  I  0  (X,x,  v,z)dxdy  (9) 

2c  2X  lc  -i-2X 

If  instead  the  averaging  is  done  only  over  the  chord  of  the  finite  wing, 
then  a  correction  will  be  obtained  which  will  be  a  function  of  aspect  ratio 
X  and  soan  y,  +c 

Bx(X»y)"?c  J  Vx»x»y»*>d*  (9) 

-c 
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With  the  correction  as  a  function  of  spanvise  position,  it  becomes 
applicable  to  a  ro4 ary  wing  where  there  are  major  spanvise  variations. 

The  variable  z,  in  the  airfoil  thickness  direction,  is  eliminated  since  the 
wing  is  assumed  to  be  in  the  z  =»  0  plane. 


THE  CORRECTION  FOR  A  SEMI-INFINITE  WING 


To  y.n-*- 


The  main  effect  of  a  complementary 
wing  will  be  felt  near  '•■he  tip 
where  the  eomp.lementarv  wing  joins 
the  finite  wing,  rather  than  at 
the  opposite  wing  tip.  7hU9  a 
further  simplification  may  be 
found  if  the  complementary  wing 
opposite  to  the  tin  under  consider¬ 
ation  is  neglected.  The  problem 
may  then  be  posed  as  the  complemen¬ 
tary  wing  of  a  semi-infinite  wing 
(see  Figure  ?.) .  Mathematically, 
this  removes  the  first  terra  of 
equation  (M,  In  this  section 
the  correction  for  the  semi- 
infinite  wing  will  be  derived, 
while  in  the  next  two  sections  the 
correction  for  a  finite  wing  is 
derived,  and  then  the  two  correc¬ 
tions  are  compared, 


Mathematically,  the  problem  is  the 
substitution  of  equation  (2)  for  A 

into  equation  (U)  for  the  complementary  wing  potential  and  then  in  turn 
into  equation  (9)  for  the  correction,  and  carrying  out  the  indicated 
mathematical  operations.  The  complementary  wing  potential  for  a  semi¬ 
infinite  wing  is,  from  equations  (2)  and  (U), 


0(x, y,z) 


_1_ 

2* 


+» 


l  L 


Fj-(0  d* 


(10) 


/  (x-O**  (y«o)1/+ 


Taking  the  derivative  of  G  with  respect  to  x  results  in 

,+C  ,  Ff;(r.)  d£  dq 


ao 

3x 


3x 


J 0  f_c  2"  Ax-O2  *  (y-R)^  *2 

»  +c 

J_  j  _  F&(0  (x-OdC  dn 

2lT  \  ■'•c  {  () 


(11) 


■O2  +  (y-n)2  + 
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The  integration  with  resrect  to  n  over  the  span  (o<n<»)  0f  the  comolenen- 
tary  wing  through  the  use  of  integral  form  Petit  Bo-fs  67-1442* 


du 

+  u2)3'2 


_ u _ 

A2  J~ A2  +  u‘ 


(12) 


where  u  ■  y  -  n 
du  “  -  dn 

A2  «  (x-O2  +  z2 

and  F  (5)  is  constant  with  respect  to  the  integration.  Noting  that  the 
upoer 'limit  of  integration  is 


u 


lim  _ 

u-*-  — »  /a2  +  u2 


-  1 


results  in 


.  +c 


x  2tt 


-c 


F5  (x-O  d£ 


+c 


aw 

r  -du 

I  (A2  +  u2) 


3/2 


x  2  T 


-1  f  r^x-adr 

(x-O2  +  z2 


J 

-c 


/(x-O2  +  y2  +  z2 


+  l 


(13) 


Substituting  into  equation  (h)  to  obtain  the  average  chordwise  correction, 


C,(v)  - 
X  4ttc 


+c  +c 

r 


TV  (x-O 


(»-«’  ♦  [  CohTT? 

-c  — c 

Evaluating  at  the  chord  line,  z  »  0,  equation  (.14)  simplifies  to 


d£  dx 
(14) 


r+c  +c 

_ 

f  r  fs  . 

1 

+  1 

1  (x-O 

/(x-O2  +  y2 

-c  -c 

“ 

d£  dx  (15) 


*  Petit  Bois  67-14  indicates  the  14th  integral  form  given  in  Reference  42 
on  page  67. 
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Equation  (15)  may  be  integrated  with  respect  to  x  holding  F^  as  a  constant, 
since  Fj-  is  a  function  of  the  auxiliary  coordinate  £.  The  use  of  integral 
forms  Petit  Bois  1-U  and  U5-10  results  in  the  following  equation: 


+c 


-1 


1  Fru) 

„„  i  +V  (^>2- 1 

1  s 

-c 

i  +V(^.)2  + 1  . 

V 

d* 


(16) 


To  evaluate  equation  (l6)  further,  a  specific  airfoil  under  consideration 
must  be  given  and  its  slope  Ft-  inserted  into  equation  (l6)  before  carrying 
out  the  integration.  An  attempt  was  made  to  carry  out  this  integration  for 
an  airfoil  whose  thickness  is  given  by  a  simple  algebraic  series,  such  as 
the  NACA  fou--digit  series,  without  success.  Two  other  approaches  were 
used  successfully  to  evaluate  the  integral  of  equation  (l6).  The  first  is 
numerical  integration,  while  the  second  is  a  Taylor  series  expansion  which 
results  in  a  simple  but  accurate  approximate  equation  of  wide  application. 

burner ical  Integration 

Numerical  integration  of  equation  (16)  va9  made  using  Simpson's  rule  in  a 
computer  program.  This,  of  course,  requires  an  even  spacing  of  intervals. 

In  most  instances,  30  intervals  were  used  with  the  upper  and  lower  limits 
+  0.9999c  to  avoid  possible  singularities  and  division  by  zero  within  the 
computer.  To  indicate  the  accuracy  with  the  number  of  intervals,  ''’able  I 
compares  the  values  of  ox  using  20  and  30  terns  for  a  parabolic  arc  airfoil. 
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Differences  in  the  final  change  in  torque  due  to  the  number  of  terms  are 
even  smaller  than  indicated  in  the  cable. 


Further  discussion  of  results  will  be  in  comparison  with  other  methods  of 
evaluating  equation  (16). 


Taylor  Series  Approximation 

If  the  bracketed  term  of  equation  (16)  is  expanded  in  a  Taylor  series  of  £, 
the  airfoil  slope  F-  will  then  appear  in  integrals  of  the  form 


FcUKnd  C 


(17a) 


which,  as  will  be  shown  later,  may  be  further  reduced  to 

+c 

F(CKm<iC  (17b) 

-c 

Anderson  found1  that  for  liie  usual  airfoil  section,  the  values  of  these 
integrals  cannot  have  large  variations. 


Therefore,  we  can  write  equation  (16)  as 
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+c 


G„=  — —  H  (c/y)  +  — H'(c/y)  I  F-ZdZ  +  ...  + - - - 

x  ;rtc  0  4;rc  1  |  ■>  4-nc  n.  | 


1  LI  ' 


H  n' (c/y) 


-be 


(18) 


-c 


-c 


where 


»«)  -  In  1  +  '  “ggfgflj-  -  Hc  +  {»•  (c/,,5  -  0) 

1  +  /TCc+O/y)2  +  1 


+  ...  4-c"  Hn'(c/y,  5-  0) 

n! 


(19) 


where 


3% 

Hn 


(20) 


The  integrals  of  equation  (18)  will  be  discussed  in  a  later  section  on  the 
effect  of  airfoil  section,  while  here  the  first  three  nonzero  terms  of  the 
Taylor  series  of  equation  (19)  will  be  given. 

For  the  Taylor  series,  the  derivatives  of  the  logarithm  terms,  which  are  of 
the  form 


tn(l  +  /~p2+  1  ) 


(2.1) 


are  found  as  follows : 
First  Derivative 


dZ 


-(in  (1  +  J  pz  • +  1 


)  = 


1 

v'p2  +  1  +  1 


1_ 

2 


2p 


p_l  i  -  + 1 

p  "ftt- 


(22a) 


J5 


Second  Derivative 


F 


£  (_  P'  y 

dC  P  /  p^+i 


*  P 


oP*/f)^  +  1  (~p/  /p2+i)-(l“  »p2+l)  (p2/ /p2+l  +  i^+l  ) 


p2(p2+i) 


’12| 


2p2+ 1 


^,-2^3/2 


-  u 


P^PM-l) 


(22b) 


Third  Derivative 


r.  (?i 

\  2P2+1  _ 1 ) 

=  -  p' 3 

p(p2+l)4p- (2p2+l)  (2(p2+l)  +  3p)  2 

i  1° 

p2(p2+l)3/5 

p4(p2+l)3  p3 

P'l  [  6p4+5p2+2  _ 


P|  l  (p2+l)5/2 


(22c) 


Fourth  Derivative 


,  f 

i  i 


d 

d<- 


=  P 


P  ’ 3  /  6p4+5p2+  2  +  2j 

p3  \  (p2+l)S;2 

p3(p2+l)/2(2/.p3+10)  - 


5/2  3/2 

(6p4+5p2+  2)  (3p2  (p2-tl)  t-  pJ  5 (p2  +  1)  p)  6 


p6(p2+l)5 


f  I 

f  .* 

i  i 


3p 


•  4 


8p6  +  8p4  +  7p2  +  2 
( p2+i)7/ * 


-  2 


(22d) 


Fifth  Derivative 


_d 

df. 


3p' 4 / 8p6  +  8p4  +  7p2  +  2  _  \ 


(p2+  1) 


7/2 
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r 


p4  (p2  +  l)  An)-(8p6+8p4-F7p2+2)  (4p3 (p?+l)  7',%-pu ?(p2-t-i)  £Ap)  8 

op  - : - - - - - — — — - +  — 

p8(p2+l)7  P5 


v5  r 

p5  L‘ 


AOCp^+pfe)  +  b  3p‘4+36p2+  8 

(p2-H)9/2 


(22e) 


The  terns  of  the  Taylor  series  are  then  found  by  evaluating  at  £  =  0.  Note 
that  there  are  two  terns  of  the  form  of  equation  (21)  with  p^  =  (c-C)/y  and 
*  (  r*-£)/y.  Further  note  that 

dp.  dp. 

K  =-dT  =  'dr  =  'p2  =-1/y 


First  Term 


H0(c/y,C  =  0) 


1  +  /(c/y)2+  1 

in  - - - —  =  in  1  =0 

1  +  /  (c/y)2  +1 


(23a) 


Second  Term 


H[(c/y,5  =  0)= 


(1/y)  1  -  /(c/y)2  +  T  ___  (-l/y)  1  -  /(c/y)2  +  1 
(c/y)  /(c/y) 2  +  1  +  (c/y)  /(c/yP  +  1 


'(c/y)2  +  1-1 
/(c/y)2  +  1  " 


(23b) 


Third  Term 


(l/y)2  2(c/y)2  4-  1  (~  l/y)2  2(c/y)'>  +  1 

y’^  (c/y)2  (  (c/y)2  +  l)3/2  (c/y)2  (  (c/y)2  +  l)3/2 


(c/y, 4  =  u;  =  0 


(23c) 


Fourth  Term 


H'"  (c/y  ,4  =  0)  =  4-2 


6(c/y)4  +  5 (c/y) 2  +  2 
(  (c/y)2  +  l)5'2 


(23d) 
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Fifth  Tent 


H,,v(c/y,C  -  0)  =  0 


Sixth  Term 


H  (c/y,£=0)  =  -5L8  -  6 


1 


40 


[(c/y)a  +  (c/y)*]  +  63  (c/y) 4  +  36(c/y)?  +  8 


(  (c/y)2  +  1  ) 


9/2 


(?3e> 


Therefore,  the  correction  factor  may  be  written  in  terms  of  a  Taylor  series 
expansion  for  the  first  three  nonzero  terms  as 


Gx  =  (I/^ttc)  (h^I^/c  +  h^^/c3  +  h3I3/c5) 
where  we  have  defined 

hx  =  H’(c/y,5  =  0) 

h2  =  H"'(c/y  ,C  -  0 ) / 3 ! 

h  =  HV(c/v,£  «  0 ) / 5 ! 

+c 


h  -/ 

-c 

-fc 

12  =J  F^3dC 

t3  ./%f55dC 


(24) 


(25) 


(26) 


'-c 


THE  CORRECTION  FOR  A  FINITE  WING 


The  derivation  for  a  finite  wing  follows  a  very  similar  pattern  to  that 
presented  for  the  semi-infinite  wing.  The  complementary  potential  for  a 
finite  wing  is 

-2A 

f  F.  df  dn 

J  k'(x-Oz  +  (v-n)^  +  z2 

+c 


G(A,x,y,z)  =  -  — 


dC  dn 

/(x-£;)^  +  (y-'l) 2  + 


(27) 
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Differentiating  with  respect  to  x, 
+c  -2X 

f 


C-  ~  -T~- 
x  2tt 


Fr  (v-r)  dr  dll 


J  (  (x-O7  +  (y-n)2  +  22)3/2 


-c 


+ 


Fr(x-0  dC  ^ 

s 

(  (x-O2  +  (y-n)2  +  z2)3^ 


(28) 


Integrating  with'  respect  to  n  over  (-“  <  n  <-2X)  and  (0<  n  <«  )  in  the  same 
manner  as  equation  (13)  was  obtained; 


-1 

GX  ~  2  71 


f  ^ 

y 

2X-y 

-  +2 

I  (x-O 

/(x-O2  +  y2 

* 

/(x-O^  +  (2X-y)2 

-c 


d( 


(29) 


Averaging  in  accordance  with  equation  (9)  and  evaluating  at  z  =  0, 
+c  +c 


Gx 

A  4nc 


r  _ 

J  (x 


•5) 


-c  -c 


2X-y 


/(x-£)2  ■»■  y"  /(x-O2  +  (2>-y)Y 


+2 


d£  (30) 


Carrying  out  the  integration  in  the  manner  used  to  obtain  equation  (16) 
results  in 

+c 


Gx.(y)  =  4nc 


1  +  /u2  +  1  1  +  Ait  +  l 

in - +  in-  - 


i  ♦  7u|~TT 


-c 


1  +■  Ai£  +  1 


d£ 


(31) 


where 


u1  -  (c-O/y 
u2  =  (c+O/y 
u3  =  (c-U/(y-2X) 
u^  =  (c+C)/(y-2X) 


S  (32) 
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We  may  compare  equation  (31)  for  the  finite  wing  with  equation  (16)  for  the 
semi-infinite  wing.  As  may  be  expected,  the  collection  factor  for  the 
finite  wing  is  the  sum  of  the  correction  factor  for  a  semi-infinite  wing 
with  its  tip  at  the  origin  going  toward  positive  infinity,  given  by  the 
first  term  (u^  and  u2)  of  equation  (31)  and  another  semi-infinite  wing 
with  its  tip  at  y  ■  -2,\,  going  toward  negative  infinity,  given  by  the 
second  term  (u^  and  u^)  of  equation  (31).  In  a  similar  manner,  the  Taylor 
series  expansion  terms  given  by  equations  (23)  are 


H 


i 


-H  (c/y 
finite  I 

wing  «SI 


c/y)  +■  (c/y  =  c/(2X-y)  ) 

SI 


(33) 


where  SI  indicates  semi-infinite  wing.  The  forms  of  equations  (24)  and  (25) 
will  chus  hold  for  finite  wings. 


AIRFOIL  SECTION  CHARACTERISTICS 


The  integrals  of  equation  (26)  may  be  simplified  and  use  made  of  the 
equation  for  the  atrfoil9  rather  than  the  airfoil  slope,  by  setting  up 
these  integrals  using  the  method  of  integration  by  parts  (/udv  -  uv-/vdu). 
The  first  integral  in  equation  (26) 

+c  +c  +c  +c 


IL  =  /  F^ds  =  SF 


Fd£ 


Fd£  (34a) 


-c 


-c  -c 


-c 


since 


1. 


dF 

d£ 


dv  =  F^dC 


dF 


du  =  d£  v  *»  F 


3.  F  •  0  at  4  •  j;  c,  i.e.,the  thickness  is  zero  at  both 
leading  and  the  trailing  edges. 


Statements  1  and  3  will  hold  exactly  for  the  other  tw^  integrals,  and 
statement  2  will  be  similar  with  u  =  ?  and  du  =  nC  df,. 


+c 


+c 


and 


(34b) 


(34c) 


The  limitation  on  these  integrals  may  be  found  from  comparing  the  actual 
profile  to  a  rectangle  enclosing  the  profile.  Comparing  the  areas  of  the 
rectangle  and  the  profile, 

+c 


2 


F  dC  =  area  of  airfoil 


-c 


Let 


then 


4tc  a  area  of  rectangle 
+  c 


6A  =  2  Fd^/4tc 


-c 


+c 


i  ,1  c  c  „  t 

-I  =  Fdf;  =— 6.4tc  —  — -  =  26,  —  c.c 
1  2  A  cc  A  c 


-c 


(35a) 


(35b) 


Nondimensionalizing  by  defining  the  leading  and  trailing  edges  at  +  1 
(chord  length  =2,  c  =  1) , 

■  “2oat  (35c) 
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Figure  3.  Geometrical  Relationship. 


The  definiton  of  6^  is  the 
ratio  of  the  area  of  the  profile 
to  the  area  of  the  rectangle.  It 
may  be  seen  from  Figure  3  that 
the  ratio  of  these  two  areas 
will  not  vary  greatly  for  the 
usual  airfoil  sections.  In  a 
similar  manner,  the  moment  of 
inertia  of  the  profile  and  the 
enclosing  rectangle  may  also  be 
compared: 


2  F  d£  “  moment  of  inertia  of  the  profile 

(2t)(2c)3/12  *  moment  of  inertia  of  rectangle 
Let 

-  jf®  2  F  5  dS  /  (4tc3/3)  (36a) 

When  nondimensionali2ed ,  Ij  becomes 

I-  *  /C  F  (3  d(  »  -  26  t/3  (36b) 

Z  C  1 


Again  by  inspection  of  Figure  3,  the  ratio  of  the  moment  of  inertia  of  the 
profile  to  that  of  the  rectangle  will  not  vary  by  much.  Unfortunately, 
there  is  no  similar  physical  argument  for  the  third  integral. 

The  values  of  the  integrals  of  equation  (34)  were  calculated  for  some 
airfoils.  Table  II  shows  these  calculated  values  for  an  airfoil  approximated 
by  an  ellipse  for  the  nose  section  and  a  parabolic  arc  for  the  rear 
portion  (from  reference  1),  and  a  parabolic  arc  airfoil. 

It  should  be  pointed  out  that  6.  and  6j  are  independent  of  the  coordinate 
system  and  are  dimensionless,  wftile  the  values  of  the  integrals  do  depend 
on  the  coordinate  system,  and  the  values  given  in  the  table  should  be  used 
only  for  the  case  where  the  length  is  measured  in  terms  of  the  semichord 
as  is  done  in  this  report  (see  Figure  3.) 
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A  simple  computer  program  was  written  to  calculate  Gx  so  that  its 
characteristics  may  be  studied.  The  program  Integrated  equations  (16)  and 
(31)  numerically  for  Cx  and  calculated  as  a  Taylor  series  expansion 
(see  equation  (24)  for  both  the  finite  and  the  semi-infinite  wing).  From 
these  calculations,  the  plots  presented  in  Figures  4  through  8  were  made. 

The  comparisons  made  in  this  section  are  for  a  10%  thick  parabolic  arc 
airfoil  section,  but  the  general  trends  will  apply  for  other  sections. 

Figures  4  and  5  show  the  coefficients  of  the  Taylor  series  for  a  semi-infinite 
wing  (Figure  4)  and  wings  of  aspect  ratio  2,3  and  4.  The  following  observa¬ 
tion  may  be  made: 


-  The  coefficients  decrease  at  a  fairly  rapid  rate  near  the  tip. 


-  The  first  term  is  dominant  by  far  over  the  full  span. 

-  The  second  coefficient  becomes  negative  about  one-half  semichord 
from  the  tip,  has  a  negative  maximum  at  about  0.9  semiehord,  and 
then  becomes  asymptotic  to  zero.  For  y>2.25c,  the  absolute  value 
of  the  second  coefficient  is  less  than  one-tenth  of  the  first 
coefficient. 

-  The  third  coefficient  is  negative  in  the  range  0.25<y/c*1.2,  with  a 
negative  maximum  at  about  y  =  0.45c.  After  becoming  positive  at 
y/c^l.2,  a  maximum  occurs  at  y/c^l.S,  but  at  this  maximum  the 
ratio  of  the  first  term  to  the  third  term  is  about  1/100.  Above 
y/c>0.1,  this  term  is  less  than  one-tenth  of  the  first  term. 

-  The  second  and  the  third  terms  are  more  nearly  equal  in  magnitude, 
and  thus  one  should  not  be  neglected  in  favor  of  the  other. 

-  At  the  wing  tip,  the  limiting  values  are: 
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hx  - *  2.0 

h2 - f  0.667 

h3  - »  0.4 

-  These  observations  hold  for  the  three  finite  aspect  ratios  and  for 

the  semi-infinite  wing.  The  similarities  of  the  shapes  of  the  curves 
and  the  magnitudes  are  obvious. 

The  correction  factor  Gx  is  shown  in  Figures  6,  7  and  8.  Figures  6  and  7 
compare  the  three  means  of  calculating  the  correction  factor: 

1.  using  the  first  term  of  the  Taylor  series  only 

2.  using  the  first  three  nonzero  terms  of  the  Taylor  series 

3.  using  a  numerical  integration  scheme  with  30  terms  from  -0.9999c 
to  0.9999c  with  Simpson's  rule. 

In  addition.  Table  III  tabulates  the  data  used  for  Figure  7  and  includes  a 
column  showing  Cx  calculated  using  two  terms  of  the  Taylor  series.  Except 
within  the  0.2  semichord  near  the  tip,  all  of  the  means  of  calculating  the 
correction  factor  are  in  good  agreement.  The  three-term  approximation  is 
in  good  agreement  with  the  numerical  integration  calculation  even  closer 
to  the  tip. 

A a  would  be  expected,  the  correction  factor  decreases  rapidly  as  the  distance 
from  the  tip  is  increased.  The  correction  factor  Gx  is  the  change  in 
velocity  (AU*)  between  the  two-dimensional  and  the  three-dimensional  flows. 

At  the  tip,  for  a  10%  thick  airfoil,  AU^/U^  is  about  3.0%,  while  one 
semi chord  inboard  it  is  reduced  to  0.15%. 

The  magnitude  of  AUM/U„  is  indicative  of  the  accuracy  required  in  calculating 
Gx.  Since  AU^/U^  is  of  the  order  of  one  percent,  a  ten-percent  error  in  5X 
results  in  only  a  one-tenth  of  one  percent  error  in 

The  effect  of  aspect  ratio  on  the  correction  factor  may  be  noted  in  Figure  7. 
Here  the  aspect  ratio  effect  is  seen  on  the  magnitude  of  the  distribution 
of  G  ,  especially  at  the  wing  center  line,  where  differences  in  Gx  for  the 
various  aspect  ratios  are  seen  to  be  important  for  the  small  aspect  ratios 
considered. 

Figure  8  comparer  the  spanwise  variation  of  the  correction  factor  near  the  tip 
for  several  aspect  ratios.  It  is  seen  that  the  differences  are  rather  small. 
The  difference  between  the  aspect  ratio  20  and  the  semi-infinite  wing  could 
not  be  shown  on  the  plot  as  it  is  a  difference  of  one  in  the  third  significant 
place.  The  effect  of  aspect  ratio  on  the  value  of  the  correction  factor 
near  the  tip  is  thus  small,  as  any  effect  is  due  to  the  complementary  wing 
associated  with  the  opposite  tip,  and  since  it  is  far  away,  it  has  only  a 
small  effect. 
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Figure  4.  Coefficients  of  the  Taylor  Series  for  a  Semi-Infinite  Wing. 


0.40 
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orrection  Factor  for  a  Semi-Infinite  Wing  With  a  Parabolic 
re  Airfoil  Section;  Thickness  =  10 Z. 


TABLE  ID.  CORRECTION  FACTOR  FOR  A  SEMI-INFINITE  WING 


SPAN 

TERMS 

FIRST 

OF  THE  TAYLOI 

SECOND 

I  SERIES 

THIRD 

INTEGRATED 
USING  30  TERMS 

"•*,  (J05M800 

-0 

,02111497 

•0,02529905 

•0,02706850 

-0,2293451* 

0,01000000 

•0 

,02100849 

"0 ,02512531 

•0,02685331 

.0,22670782 

0.01500000 

•  0 

,02090240 

•0,02499562 

•0,02663823 

, 2261808? 

0x02000000 

-0 

,02079839 

•0 , 02478600 

•0,02642328 

.-‘,22770379 

0,02900000 

•0 

,02069033 

•0,02461648 

-0, 02620853 

-.’,22725842 

0,09000000 

-0 

,02016097 

•0,02377112 

■0,02513926 

-'’,.’2531807 

0,10000000 

*0 

,01910914 

•0 , 02210093 

•0.02304137 

-7,72226477 

0,19000000 

-0 

,01807279 

•0,02047623 

•0,02103217 

-7  ,31986488 

0,20000000 

-0 

,01709896 

-0 ,01891626 

•0,01914564 

-7, .’1788802 

0,29000000 

•  0 

,01607391 

•0,01743656 

■0,01740549 

-7, .’1621633 

0,30000000 

,01512296 

•0,01604843 

■0, ? 1582454 

-.’.,‘1477810 

0,39000000 

,01421042 

•0,01475888 

•0,01440564 

-7,01352551 

0,40000000 

•0 

,01333952 

•0,01387091 

•P, 01314358 

-7,31242437 

0,49000000 

-0 

,01251246 

•0,01248410 

•0.01202747 

•7,01144918 

0,90000000 

-0 

,01173050 

•0,01149527 

•0,01104306 
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Thus  a  result  to  be  found  later  may  be  anticipated.  In  a  rotor,  which  is 
a  blade  of  large  (on  the  order  of  15  to  20)  aspect  ratio,  a  large  part  of  the 
forces  are  near  the  tip,  and  this  is  where  roost  of  the  correction  needs  to 
be  applied.  But  near  the  tip,  the  actual  aspect  ratio  has  a  second-order 
effect.  Therefore,  for  rotary-wing  calculations,  the  simplified  semi- 
infinite  wing  correction  factor  will  be  found  satisfactory. 

CORRECTION  FACTOR  FOR  A  ROTATING  WING 

Sears,12  in  1950,  derived  the  relationship  between  the  incompressible  veloci¬ 
ties  on  a  rotating  blade  and  on  a  similar  cylinder  in  plane  steady  flow. 

The  relation  of  the  potentials  is 

♦  -  fty^  -  x)  (37) 

where  4,  is  the  potential  of  the  rotating  blade  and  ^  is  the  potential  of 
the  blade  in  uniform  flow  at  unit  velocity. 


The  velocities  in  the  transverse  (to  the  span)  plane  are  identical  in  the  two 
cases.  The  spanwise  component  of  the  velocity  is  given  by  a  very  simple 
formula,  and  for  an  infinite  circular  cylinder  this  velocity  is  zero  on  the 
surface  and  therefore  makes  no  contribution  to  the  forces. 

A  recent  paper  extending  Sears  by  Goorjian  and  McCroskey13  indicates  that 
for  finite  span  rotating  blades,  there  is  a  finite  spanwise  potential 
velocity  at  the  surface.  While  this  component  is  important  in  boundary  layer 
calculations  and  thus  affects  the  separation  line,  its  contribution  to  the 
pressure  in  the  linearized  case  being  applied  here  is  of  the  same  order  of 
magnitude  as  other  terms  which  are  neglected. 


Therefore,  for  a  rotating  blade,  we  may  write, 
discussion  for  fixed  wings. 


u  =  - 


fiy  + 


/ 


plan form 


fly 

2ir 


from  equation  (37)  and  the 


(x-g)dS 

F5  d3 


(38) 


where 


D2  -  (x-5)2  +  (y-n)2  +  z2 


(39) 


and 


dS  -  d(dri 


(40) 


In  a  manner  similar  to  that  used  for  the  fixed-wing  case,  the  integration 
over  the  finite  planform  may  be  split  up  as  integration  over  the  infinite 
blade  less  the  integral  over  the  complementary  blade.  Averaging  over  the 
chord  results  in 
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where  ^  is  the  velocicy  due  to  the  infinite  blade.  The  spanwise  coordi¬ 

nate  y  in-equation  (38)  may  be  brought  outside  of  the  integral  sign  since  it 
is  invariant  with  respect  to  each  of  the  three  variables  of  integrations. 

This  equation  is  identical  to  the  result  obtained  for  the  serai-infinite  wing 
(equation  (6b)  with  the  averaging  given  by  equation  (9)  and  the  integrand  in 
equation  (11)).  Thus  it  has  been  shown  that  the  correction  factor  in  the 
case  of  a  rotating  blade  is  given  by  the  same  equation  as  for  a  fixed  wing. 


CORRECTION  FO  THE  CHANCE  IN  DYNAMIC  PRESSURE 


The  correction  factor  is  the  difference  between  the  velocities  of  a  finite 
wing  and  an  infinite  wing*  having  the  same  potential.  This  in  turn  means 
that  the  two  wings  have  the  same  static- to-total-pressure  ratio  and  therefore 
the  same  force  acting, 


finite  wing  infinite  wing 

with  the  velocity  condition  of  equation  (7),  which  requires  different 
dynamic  pressures.  Thus,  since  the  reference  areas  are  also  the  same. 


(A2) 


(CDq)FW  =  (CDq)A  = 
but  neither  CQ  nor  q  is  the  same.  Thus 


(43) 


(44) 


or  by  logarithmic  differentiation  of  (43) 

^CD/CD  =  “  A  q/q  (45) 

The  correction  gives  the  change  in  velocity  and  thus  the  change  in  dynamic 
pressure.  Since  the  drag  is  a  function  of  Mach  number,  the  change  in  velocity 
must  be  applied  not  only  as  a  change  in  drag  coefficient  but  also  as  a  change 
in  Mach  number  as  will  be  shown  in  Figure  10.  It  is  more  convenient  to 
express  both  of  these  changes  in  terms  of  the  f rcc-streara  Mach  number. 


In  the  introduction,  the  explanation  of  the  tip  relief  effect  was  made  on 
the  basis  of  deviation  allowed  of  the  stream.  Consider  a  representative 
stream  tube  as  it  passes  over  a  finite  and  an  infinite  wing  illustrated  in 

*  While  the  discussion  here  talks  about  wings,  since  the  correction  factor 
was  derived  as  it  varies  along  the  span  of  a  finite  wing,  the  discussion 
presented  here  is  just  as  valid  for  a  spanwise  element  of  the  wing  located 
a  given  distance  y  from  the  tip. 


32 


Figure  9.  Distortion  of  a  Stream  Tube  Over  a  Finite  and  an  Infinite  Wing. 

Figure  9.  In  the  case  of  the  infinite  wing,  the  stream  tube  cross-sectional 
area  is  decreased  and  displaced  due  to  the  displacement  effect  of  the  wing. 
But  for  the  finite  wing,  tee  stream  tube  may  also  expand  laterally;  thus  the 
cross-sectional  area  is  decreased  by  a  lesser  amount. 


The  flow  in  both  of  these  stream  tubes  is  isentropic,  and  thus  the  usual 
quasi-one-dimeusional  isentropic,  compressible  flow  analysis  may  be  used. 
Since  both  stream  tubes  have  identical  conditions  at  infinity,  they  may  be 
considered  as  part  of  the  same  stream  tube.  In  the  usual  application  of 
this  analysis,  the  area  changes  are  known.  In  this  case,  however,  the 
velocity  changes  are  given  by  the  correction  factor,  while  the  area  change 
in  quantitative  terms  is  immaterial.  It  is  desired  to  find  AC^/Cy  and 
i.M/M  in  terms  of  £U/U  =-Cv,  with  the  appropriate  compressible  transformation. 

The  equation  of  continuity  for  a  stream  tube  may  be  written  as  either 

PUA  =  Const  (46) 

or 
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where  A  is  the  area  and  the  starred  quantities  represent  conditions  at  a 
sonic  point  as  a  reference.  Also  Bernoulli's  equation  is 


lid  Li  +  dp/p  =  0 

Logarithmic  differentiation  of  equation  (46) 


3p  3U  3A 

—  +  ~T  +  ~T  =  0 


combined  with  Bernoulli's  equation  leads  to  the  familiar 


(M2-  D-SL.  dA 

u  1  U  A 


(48) 


(49) 


(50) 


used  so  often  in  describing  the  difference  between  subsonic  and  supersonic 
flow.  Using  the  relations  from  iser.cropic  flow 
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Equation  (47)  results  in,  after  some  manipulation, 
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=  constant 
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Logarithmic  differentiation  gives 
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(52) 


(53) 
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Combine  equations  (49)  and  (53) : 
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dM  Y  -  1  , 

— •  (1  +— «2> 


dU 

U 


(54) 


which  is  the  relation  between  the  velocity  change  as  given  by  the  correction 
factor  and  the  resulting  Mach  number  change. 

Since  q  «»  pU 2/2,  logarithmic  differentiation  gives 


dq  dp  dU 

q  ~  p  +  2  u 


(55) 


Combining  with  equation  (49)  gives 


dq 
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dA  dU  dU  dU  dA 
A  “  U  +  2U  =  U  ”  A 


(56) 


Substituting  from  equation  (50)  and  then  from  equations  (54)  and  (45), 
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(57) 
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(58) 


As  previously  noted,  the  derivation  here  is  applicable  to  a  segment  of  a 
wing  or  to  a  full  wing,  depending  on  how  the  correction  factor  is  derived. 

In  fact,  equations  (54)  and  (58)  are  given  by  Anderson.1  But  since  he  did 
not  give  their  derivation  in  detail  or  explain  them  in  a  physical  manner, 
and  since  they  are  unusual  in  form,  it  is  felt  that  presenting  their 
derivation  in  detail  would  be  useful  in  clearing  up  some  misunderstandings. 

Figure  10  illustrates  the  application  of  equations  (54)  and  (58).  The 
correction  factor  used  here  was  from  Reference  1  and  for  aspect  ratio  3  wing 
with  NACA  0012  airfoil  with  two-dimensional  drag  vs.  Mach  number 
characteristics  approximated  by  two  straight  lines  as  proposed  in  Reference  55. 
A  point  on  the  two-dimensional  curve  is  displaced  horizontally  in  accordance 
with  equation  (54)  and  then  vertically  as  per  equation  (58)  to  obtain  the 
three-dimensional  characteristic. 

In  reading  this  plot,  the  drag  of  the  three-dimensional  wing  is  seen  to  be 
less  than  the  drag  of  the  corresponding  two-dimensional  wing  at  the  same 
Mach  number. 
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Experimental  verification  of  the  basic  theory  for  fixed  wings  may  be  found 
in  the  report  by  Anderson  and  Carroll.43  Wings  of  various  aspect  ratios 
at  high  subsonic  speeds  are  correlated  both  for  the  case  of  nonlifting 
wing  and  for  the  lifting  case  with  lift  coefficients  up  to  0.3. 

The  corrected  drag  coefficient  as  calculated  for  Figure  10,  requires  the 
use  of  appropriate  transformation  between  the  incompressible  theory 
as  derived  here  and  the  compressible  flow  considered.  For  this  fixed  wing 
case  the  direct  Prandtl-Glauert  transformation  is  used.  For  rotary  wing, 
the  extension  discussed  in  the  next  chapter  must  be  made  to  apply  the 
transformation  shown  in  equations  (105)  to  (110). 


Figure  10.  Theoretical  Correction  of  the  Drag  Coefficient  vs.  Mach  Number 
for  an  HACA  0012  Airfoil  From  Two-Dimensional  Data  to  an 
Aspect  Ratio  of  3.  (Note  Displaced  Origin.) 
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BASIC  EQUATIONS  OF  A  ROTOR  IN  COMPRESSIBLE  FLOW 


Having  derived  a  tip  relief  correction  for  a  fixed  wing,  the  next  step  is 
the  problem  of  applying  this  correction  to  a  helicopter  rotor.  To  do  this, 
the  basic  fluid  mechanic  equations  of  a  rotor  have  to  be  examined.  There¬ 
fore,  in  this  section  the  Navicr-Stok.es  equations  for  a  rotor  in  compressible 
flow  will  be  written. 

The  usual  approach  to  compressible  flows  is  the  small  perturbation  theory, 
and  this  derivation  will  be  presented.  In  fixed-wing  usage,  the  result  is 
the  Prandtl-Glauert  transformation.  It  will  be  shown  here  that  further 
assumptions  are  needed  so  that  a  Prandtl-Glauert  type  transformation  between 
compressible  and  incompressible  flows  may  be  made  for  a  rotor. 

Previous  work  on  the  representation  of  a  rotor  in  compressible  flow  by  a 
system  of  sources  and  sinks  was  done  by  Sopher.3  This  work  will  be 
thoroughly  reviewed  to  show  how  it  is  connected  to  the  present  presentation. 
The  similarities  and  differences  will  be  discussed,  along  with  where  each 
approach  is  useful. 

FLUID  DYNAMIC  EQUATIONS  FOR  A  ROTOR 

Consider  a  coordinate  systeu  with  a  fixed  origin  0  rotating  with  angular 
velocity  Q.  The  absolute  velocity  is  while  the  velocity  relative  to  the 
rotating  frame  is  tf'.*  The  relation  between  the  two  velocities  is 

V  «  V‘v  n  x  r  (59) 

The  equations  relative  to  the  rotating  frame  are  the  same  as  for  a  fixed 
frame  except  for  the  acceleration  term 


3t 


+  v ' . vv  +  a  x  v 


With  some  vector  identities,  we  have 

(V'*V)  V  =  (V’-V)  V'  +  H  x  V' 
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■—  x  r)  =  0 
3r 


(60) 
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*  This  discussion  is  taken,  to  a  large  extent,  from  Hughes  and  Gaylord.44 
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Therefore,  Che  equation  of  momentum  is 


— -+  (V  V)V'  +  2«  x  V’+  5?  x  (ft  xr)  -  -  Vp  (62) 

where  the  last  two  terms  on  the  left-hand  side  are  due  to  the  Coriolis  and 
the  centrifugal  forces.  Let  the  velocity  have  components  vf,  vet,v  in 
a  cylindrical  reference  frame,  and  define  the  axis  such  that  rotation  fs 
only  about  the  z-axis.  Further  assuming,  for  convenience,  steady  flow 
with  rotating  coordinates,  the  three  components  of  the  momentum  equation 
are 


va  3v 

r  30 

3vr 
Vz  3z 

v0t 

~  -  2V  ‘ 

rQ2 

P 

i£ 

3r 

(63a) 

3v0t  V0t 
—  - 

3v0t 

3vr  veyr 

2vrft  -  ~j-r 

3p 

(63b) 

3r  r 

30  +  V 

2  “r  +  r“  + 

30 

3vz 

\T  ' 

vA.  3v„ 

0t  z 

3vz  1 

ap 

r  3r 

-r  t  v  - 

r  39  z 

3z  p 

3z 

wJCj 

The  equation  of  continuity  is  identical  in  rotating  and  stationary  frames 
of  reference: 


3v  i  3v 

r  _1_  et 

3r  r  36 


3v_  v  v_ 

+  2  +  _L  +  — £ 

3  7.  r  p 


v0t3p  Vz  3p 
3r  +  Pr  39  +  p  3z 


(64) 


Further,  the  speed -of-sound  relation  is 

(3p/3p)  *  ; 


(65) 


To  combine  equations  (63)  ,  (64)  and  (65)  by  multiplying  each  of  the 
equations  in  (63)  by  v  ,  vytand  v  ,  and  by  using  equations  (64)  and  (65) 
to  put  the  pressure  an5  the  density  in  terms  of  the  speed  of  sound  gives 
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Because  of  the  choice  of  axis,  in  terms  of  the  its  compoients  V  can  be  written 
as 


V  =  vr6r  +  (vQ  -  nr)§Q  +  vz6z  (67) 

where  er,  eQ ,  e2  are  unit  vectors.  The  theta  component  is  made  up  of  two 
parts,  one  a  "free  stream"  of  magnitude  -Qr  and  a  difference  velocity  v0, 
which  together  make  up  v0c.  The  only  rotational  component  in  equation  (9) 
is  the  nr  component,  thus  we  may  write 
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With  these  two  conditions,  equation  (66)  may  be  written  as 
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The  energy  equation 


a2  +  ((y  -  l)/2)  V  •  V  =  Constant  (72) 

may  be  used  to  relate  the  velocity  components  and  the  reference  speed  of 
sound  a»  to  the  local  speed  of  sound  a.  The  reference  speed  of  sound  is 
the  speed  of  sound  in  the  undisturbed  regie"  of  the  fixed  frame.  In  the 
rotating  frame,  this  is  the  speed  of  sound  at  velocity  -Cr;  thus 

a2  +I^i(Qr)2  =  a2  +  v’  •  V’ 
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If  we  also  say  that  a  velocity  potential  exists  (due  to  equation  (67)), 
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after  some  rearrangement,  equation  (71)  becomes 

$rr  +  U-  "  (Hr/a^)2}  <J>Q0/r2  +  4>zz  +  <?r/r  = 

I  (p  \ 

4>rr/a2  {-(yl)<|>0fl  +  Ip  (—•  +  <f>2  J  +  Qp 

+  $0e/ra2  {-(y+l)<f>0O  +  ~  (<Pe/r)z  +~  ^2+<i>2J} 

+  tzz^  {-(y-Dtgfi  +  ^  4r  +  Ue/r)2  +  ^  *z  ) 

Ue/r)^r  _  2 $e  »gfl  ^  2^.  /♦0£ 

a£r  ^r  a2  \  r  r2/\  r  / 


(4>e/r)  2<f.r 

M 

air 

„  Mr 

2  -tJL  <p 

+  2  — - 

_2  rz 

a* 

a* 

00 

40 


i 


It  should  be  emphasized  chat  this  equation  is  still  an  exact  equation. 

Its  counterpart  in  uniform  flow,  for  example,  would  be  equation  32-17  of 
Heaslet  and  Lomax. 14  ::1  All  terms  on  Che  right-hand  side  should  vanish  for 
the  small  perturbation  equation.  This  Is  accomplished  by  assuming  that 
the  ratio  of  all  of  the  velocity  components  to  the  speed  of  sound  is  much 
less  than  one.  Further,  it  is  assumed  that  gradients  of  the  velocities 
are  also  small,  so  that  cross  products  of  velocities  and  vejocity  gradients 
are  also  small.  Terras  of  the  rorm 

Hr  v. /a2  (77) 

I  00 

may  also  be  neglected.  The  result  is 

$rr  +  $zz  +  C1  -  ^r/aoo)2)<^QQ/rz  +  4>r/r  =  0  (78) 

In  the  case  where  ur/ao  is  nearly  one,  the  coefficient  of  the  terra 
in  the  left-hand  side  of  equation  (76)  becomes  very  small,  and  thus  the 
term  on  the  right-hand  side  must  also  be  considered.  It  is  obvious  and 
analogous  to  the  uniform  flow  case  that  the  term  of  the  form  of  equation 
(77)  is  the  largest  of  the  terms  and  we  should  write 

<j>  +  p  +  (1  -  <Qr/a  )2)  <J>aft/r2  +  $  /r  = 

rx  2z  <»  r 


-  (y  +  1)  (n/ra2M)  *e  *ee  (79) 

j 

i  For  the.  pressure,  we  have 

|  dp  =_pd(V-V/2)  (80) 

I 

Integrating  between  a  point  where  p  =  and  V  =-Qr  and  a  point  on  the 
b lade , 

^  +  * 

p  “  Poo  =^2-2~—  (  (ttr)2-  (  (ve  -  iir)2+vz2  +  Vf2  )  ) 

f 

t  Vg2+V2+V2 

I  -  -  <p-pro)(-s?rv0  + - 5 - 1  )  (81) 

i 

Again  assuming  that  all  velocities  are  small  and  that  p'<<  pro  results  in 

P  "  P~  =-pTC  ve  fir  (82a) 

or,  in  terms  of  the  pressure  coefficient. 
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The  basic  equations  (78)  and  (82) 
may  be  written  in  Cartesian 
coordinates,  as  suggested  by 
Sears12  and  defined  in  Figure  11. 
Using  che  transformations 

r2  =  x2  +  y2 

0  =  tan-1(-x/y) 

(83) 

x  =  -  r  sine 


y  -  +  r  cotO 


Figure  11.  Coordinate  System. 
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2  (yo  -  x(?  ) 

C  - - - - 2_ 

P  :;(x2  +  y2) 


(85) 


The  results  of  equation  (78)  and  (82)  are  identical  to  the  results  obtained 
by  Sopher^  directly  from  the  acoustic  equations.  Similar  work  has  been  done 
for  a  rotating  and  advancing  blade.  Simonov  and  Christianovitch'*^  derived 
the  equivalent  of  equation  (78)  in  a  similar  manner.  Also,  Bums47  derived 
the  equivalent  of  equation  (84)  directly.  We  have  gene  into  details  here  so 
that  the  further  approximations  which  need  to  be  made  may  be  better  under¬ 
stood  . 
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Using  the  transformation  6; 
acoustic  equation 


. t  in  equation  (78)  and  (82)  results  in  the 


i2  V2d>  = 


(86a) 


DIRECT  SOLUTION 


P~  P.  =  -Pa 


(86b) 


Source  in  a  Rotating  Flow 

Recently,  Sopher3  presented  the  solution  of  the  small  perturbation 
equations  using  a  source-sink  distribution  approach.  The  potential  due  to 
a  point  source  travelling  on  a  rotational  path  is 

<p  =  -k/4?rD  (87) 

where  k  is  the  strength  of  the  source  and  D  the  distance  function  referred 
to  coordinates  rotating  with  the  blade,  which  is  given  in  polar  coordinates 
by 


D2  =  rr'(0  -  6')2  +  (1  - 


)  ((r-  r')2  +  z2)  (88) 


where  the  primed  coordinates  refer  to  the  location  of  the  source  and  the 
unprimed  coordinates  are  the  field  points  where  $  is  to  be  evaluated.  In 
deriving  this  distance  function,  higher  order  terms  in  the  expansion  for 
the  cosine  are  neglected  early  in  the  derivation.  A  comparison  with  the 
more  familiar  distance  function  for  uniform  compressible  flow  in  Cartesian 
coordinates , 

Djp  =  (x  -  x’)2  +  (1  -  (y-y')2  +  z2}  (89) 

oc 

reveals  certain  similarities  and  differences.  The  form,  of  course,  is  very 
similar.  But  the  compressibility  term,  namely,  Mach  number  (U^/a^)  in  the 
uniform  flow  case,  changes  to  the  local  Mach  number  at  the  root  mean  square 
radius 


of  the  source  and  the  field  points.  Further  transformation  to  Cartesian 
coordinates  of  the  first  term  results  in  a  very  complex  expression. 

The  resemblance  between  the  two  equations  may  be  made  closer  if  it  is 
assumed  that  the  field  points  and  the  souice  points  are  close  to  each  other. 
A  further  restriction  is  that  the  two  points  are  near  the  y  axis,  which 
physically  means  narrow  blades  with  the  axis  along  the  span  of  the  blade. 
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There  are  several  arguments  as  to  why  this  should  hold,  the  most  important 
being  that,  if  these  points  are  far  apart,  D  in  equation  (87)  is  large, 
making  the  contribution  to  $  vary  small.  With  this  assumption,  it  can  be 
said  that  the  root  mean  square  radius,  r,  of  equation  (90)  does  not  differ 
greatly  from  r  or  r* .  This  will  also  provide  for  the  small-angle  expansion 
of  the  transformations  (see  equation  (83)  ),  namely. 


x  =  r  sin6  £  ^6 
v  =  -r  cose  r 


(91) 


With  this,  the  approximation  of  the  distance  function  for  rotating  flow  in 
Cartesian  coordinates  is 


D2 

Rot . 
Approx. 


=  r2  (9  -  e')2 


=  (x  -  x')2  + 


Hf)  ]c(y-y,)2 


2+ 


+ 


z2) 


z2) 


(92) 


The  resemblance  between  the  distance  functions  in  rotating  and  uniform  flow 
is  now  much  closer.  The  compressibility  factor  in  rotating  flow  is  not 
constant,  and  thus  a  "local  Mach  number"  has  to  be  defined: 

My  -  n  y/a^  (93) 

which  is  the  Mach  number  now  used  in  blade  element  rotor  analysis. 

This  derivation  of  the  approximate  distance  function  will  be  used  as  part 
of  the  justification  of  the  approximations  made  here. 

The  potential  of  a  rotor  blade  is  the  potential  due  to  the  sum  of  all  the 
sources  representing  the  blade.  This  is  the  surface  integral  over  the 
planform  area  of  the  blade. 

As  was  shown  in  equations  (82),  the  pressure  distribution  and  therefore 
the  forces  are  related  in  compressible  perturbation  flow  to  the  tangential 
velocity  vQ .  Further,  Sopher  (unpublished  notes  for  Reference  3)  derived 
from  the  usual  flow  tangency  condition  that  the  required  source  distribution 
to  describe  the  thickness  of  a  rotating  wing  is 

X  (r',8\  z  =  0)  =  -2  fir  F0/r  =  -2fi  F0  (94) 

where  F  defines  the  vertical  z  dimension  of  the  blade.  Equation  (94) 
resembles  the  linearized  uniform  flow  source  distribution,  namely,  twice 
the  free-stream  velocity  times  the  local  slope.  It  is  important  to  note  that 
the  slope  required  in  equation  (94)  is  the  slope  in  the  direction  tangential 
to  the  rotation  rather  than  normal  to  span. 
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The  integral  to  be  evaluated  then  becomes 
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«r' (6-0' )Fq ids' 


0  2tt 


D“ 


C95) 


The  surface  of  the  integration  is  the  planform  area  of  the  blade  and  in 
cylindrical  coordinates  dS'  =  r'dQ'dr'.  This  integral  was  numerically 
integrated  in  Reference  3  with  important  refinements  in  the  areas  of  using 
a  uniformly  valid  surface  flow  speed  and  an  equivalent  linear  theory  tip 
Mach  number. 


The  Complementary  Wing  Directly  in  Rota.in^  Coordinates 


An  attempt  to  apply  the  reasoning  of  the  previous  section,  the  complementary 
wing  approach  in  the  case  of  a  rotating  blade  represented  by  the  source 
distribution  in  compressible  flow,  will  now  be  made.  Because  of  the 
variable  dynamic  pressure  reference  for  the  pressure  coefficient,  the  total 
forces  are  related  to  the  total  local  velocity,  which  is  the  sum  of  fir  plus 
v.  and  will  be  denoted  by  V.  Consider  the  integration  required  in  equation 
(95);  as  previously,  split  the  integration;  from  0  to  infinity  minus  from 
tip  to  infinity.  Thus  we  can  write  the  equation  for  total  local  velocity  as 


V  -  fir  +  V0  =  fir 


fir'  (6-0')  F  ,  r'  d9'  dr'/D3 
0' 


fir'  (9-0')  F  ,  r'dO’  dr'/D3 
0 


(96) 


c  R 


where  c  represents  the  integration  over  the  constant  chord. 

In  the  case  of  the  fixed  wing,  Anderson1  obtained  a  transformation  between 
the  finite  wing  and  infinite  wings.  The  transformation  could  be  simplified 
because  of  the  following  two  results: 

(1)  The  complementary  wing  integral  (l.e.,  the  second  integration  in 
equation  (96)) reduces  to  a  form  such  that  it  is  a  multiple  of  the 
frec-stream  term. 

(2)  'The  infinite  span  integral  is  a  known  function  obtained  from  other 
theoretical  work,  from  experiments,  or  from  elsewhere. 

The  first  criterion  is  easily  met,  as  may  be  seen  by  inspection.  The  two 
terms,  after  the  integral  is  carried  out,  will  reduce  to  a  form  fi(r  +  G) 
which  may  be  used  instead  of  fir.  Thus,  we  may  write 

V  -  (fir)  effective  +  Vq  infinite  blade  ^ 
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The  problem  occurs  in  relating  the  infinite  span  integral  to  a  known 
function.  In  the  fixed-wing  problem  this  is  the  much-studied  case  of 
two-dimensional  flow  over  an  airfoil.  Many  two-dimensional,  compressible 
flow  airfoil  data  are  available.  Further,  the  Prandtl-Clauert 
transformation  is  available,  so  tnat  the  flow  may  be  compared  to  incompressi¬ 
ble  flows.  In  this  case,  though,  the  flow  over  a  rotor  of  infinite  span 
in  the  compressible  regime  must  he  known,  calculated,  or  measured  in  some 
manner.  There  appears  tc  be  no  straightforward  way  of  determining  this 
value.  Analytically,  it  brings  problems  such  as  supersonic  flows,  etc. 
Experimental  testing  is  conceptually  impossible. 

The  reason  that  this  problem  may  not  be  approached  in  the  direct  manner  is 
that  there  is  a  basic  difference  between  fixed-wing  and  rocating-wing 
flows  in  the  compressible  regime.  This  was  indicated  by  Sopher's  results3 
as  discussed  previously. 

In  incompressible  flow,  Sears 12showec  that  there  was  a  straightforward 
relation  between  uniform  flow  and  flow  over  the  same  body  rotating.  This 
same  type  of  relation  does  not  hold  for  the  derivation  shown  above  for 
compressible,  rotating  flow,  to  a  large  extent  due  to  the  fact  that  we  do 
not  yet  have  a  transformation  between  compressible  and  incompressible  flows. 
Therefore,  how  such  a  transformation  might  be  obtained  is  described  in  the 
next  section. 


PRANDTL-GLAUERT  TYPE  TRANSFORMATION  FOR  A  ROTATING  BIADE 


The  basic  differential  equations  in  Cartesian  coordinates  (equations  (84) 
and  (85))  may  be  rewritten  as 
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Consider  the  differential  equations  at  points  far  from  the  center  of  rotation. 
Now  y  is  the  spanwise  distance  and  is  of  the  order  cf  magnitude  of  the  radius. 
At  the  same  time,  x  is  of  the  order  of  magnitude  of  the  chord.  To  compare 
the  magnitude  of  each  of  the  terms  In  the  right-hand  side  of  equation  (98), 
nondimensionalize  the  x-dimension  with  the  chord  c  and  the  y-dimension  with 
radius  R.  (The  primed  coordinates  in  this  section  are  nondimens ional . ) 
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For  helLcopter  rotor  the  ratio  c/R  is  much  smaller  than  one.  The  coeffi¬ 
cients  of  all  terms  but  one,  (c/R),  are  raised  to  a  power  which  will  make 
all  terms  small  in  comparison  with  the  term  of  order  (1).  Therefore,  the 
differential  equation  becomes,  when  neglecting  the  higher  order  terms. 


(1  --  (fly / a  )^)<>  +  $  +  <p 

00  xx  yy  22 


=  0 


(101) 


In  a  similar  manner,  for  the  pressure  coefficient, 

1  c  x'  1 

2  cV  ~  ~R~  y '  ~R  $y' 


..  c2  x' 2 
"y  1  R2  ~2+  1 


(102) 


As  before,  there  are  terms  with  a  coefficient  of  one  to  compare  with  terms 
with  coefficients  (c/R)2,  which  are  small  and  neglected,  so  that 


cp  »  2<?x/fly 


(103) 


Now  equation  (101)  resembles  very  closely  the  Prandtl-Glauert  equations  of 
uniform  flow,  but  with  a  variable  Prandtl-Glauert  factor 

&2y  -  1  -  (fly/aj2  =  1  -  M2  (104) 

This  indicates  that  a  Prandtl-Glauert  type  transformation  may  be  made 
locally  for  a  rotating  blade.  The  usual  illustration  for  the  Prandtl- 
Glauert  transformation  is  by  the  change  in  chord  length  (for  example  see 
Shapiro,39  page  320).  For  a  rotor,  this  change  in  chord  length  varies 
along  the  span  due  to  the  variable  local  Mach  number.  The  resulting 
equivalent  incompressible  flow  rotor  is  illustrated  in  Figure  12. 

Since  we  now  have  a  transformation  betv;een  the  compressible  and  the 
incompressible  flows  over  a  rotor,  the  complem  ntary  wing  effect  as  derived 
in  the  previous  section  may  be  applied. 
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A  short  discussion  of  the  implications  of  this  approximation  is  in  order. 
This  approximation  is  consistent  with  blade  element  analysis,  where  an 
element  o,'  the  hlade  is  assumed  to  act  as  a  similar  airfoil  in  uniform 
two- dimensional  flow.  The  transformation  has  been  applied  previously, 
without  proof,  for  example,  by  Jones  and  Rao.45 

The  approximation  is  made  of  the  initial  equation  used  by  Sopher.3  Thus 
the  differences  found  by  Sopher  between  three-dimensional  rotating  flow 
theory  and  blade  element  analysis  in  the  local  induced  velocities  apply 
when  using  the  proposed  transformation.  Tne  problem  with  the  analysis 
presented  by  Sopher  is  one  of  application.  T^e  major  use  of  the  Sopher 
analysis  would  be  in  boundary  layer  calculations  where  it  is  important 
to  know  the  local  velocity  gradients,  and  where  Sopher 's  analysis  does 
give  a  previously  unknown  spanwise  pressure  gra.'.ien 

For  the  purposes  of  preliminary  design  and  gross  performance  prediction, 
the  approximations  made  here  are  consistent,  and  the  Prandtl-Glauert  type 
transformation  may  be  safely  used. 

EXTENSIONS  OF  THE  COMPRESSIBLE  FLOW  TRANSFORMATIONS 


The  Prandtl-Glauert  transformation  underpredicts  the  increase  in  flow  speed 
with  Mach  number  as  may  be  seen  in  Figure  C8a  of  Reference  49.  Also,  the 
Goethert  extension  to  three-dimensional  flows  considers  only  the  effect  of 
compressibility  on  lift  and  the  lift  curve  slope.  This  is  done  through 
the  change  in  wing  aspect  ratio,  since  the  Prandtl-Glauert  transformation 
changes  the  chord  but  does  not  alter  the  spanwise  or  the  thickness 
coordinates.  No  mention  is  made  of  any  effect  of  compressibility  on  the 
three-dimensional  drag. 


There  have  been  a  number  of  attempts  to  obtain  similarity  rules  that  would 
extend  the  Prandtl-Glauert  transformations  to  higher  Mach  numbers.  The 
most  important  of  these  for  two-dimensional  flow  were  reviewed  and  compared 
by  Bennett.50  Van  Dyke's  second-order  theory51 >52 »52  is  considered  to  give 
the  best  compressible  pressure  distributions  of  these  techniques,  but  its 
use  is  rather  complex.  This  method  takes  into  account  squares  and  products 
of  the  thickness,  camber  and  angle-of-attack  terms. 

But  a  method  suggested  by  Kuchemann  and  Weber54  appears  to  be  much  simpler 
to  apply.  The  usual  form  for  the  transformation  is  modified  to  the  form 

s2  -  1  -  M2(1  -  cpi) 


where  C  .  is  the  incompressible  flow  pressure  distribution.  This  form  was 
suggested  for  use  in  calculating  local  velocities  and  pressures  as  C  .  is 
a  function  of  position. 


A  simple  use  of  the  latter  transformation  was  suggested  by  Sopher'  for 
rotors;  namely,  that  an  average  value  of  Cp^  be  used,  which  can  be 
considered  constant  along  both  the  span  and  the  chord.  Thus  for  a  rotor, 
Oir/aJ/l-C  t  is  an  equivalent  Mach  number  which  may  be  used  in  the  linear 
theory. 

49 


APPLICATION  OF  THE  TIP  RELIEF  EFFECT 
IN  ROTOR  PERFORMANCE  CALCULATIONS 


Having  derived  a  correction  factor  for  the  effective  free  stream  to  relate 
the  two-dimensional  airfoil  data  to  three-dimensional  wing  characteristics, 
and  a  transformation  between  ccmpressible  and  incompressible  flow  over  a 
rotor,  we  are  now  ready  to  use  this  theory  in  helicopter  perf ormatice  cal¬ 
culations,  which  are  based  on  blade  element  analysis. 

Blade  element  analysis  involves  a  long  and  complicated  computer  program. 
Those  programs  are  generally  available  to  those  who  would  use  the  theory, 
although  it  was  not  available  for  use  within  the  present  study.  Therefore, 
the  theory  is  not  used  in  a  full  blade  element  analysis.  Presented  herein 
is  a  detailed  description  of  how  one  should  modify  existing  blade  element 
computer  programs  to  use  the  tip  relief  theory,  and  a  simplified  blade  ele¬ 
ment  analysis  for  u  hovering  rotor  using  the  theory,  so  that  an  indication 
of  the  tip  relief  effect  may  be  found. 

THE  PHYSICAL  CONCEPT  OF  TIP  RELIEF  APPLIED  TO  A  ROTOR  BLADE  ELEMENT  ANALYSIS 


The  application  of  the  tip  relief  concept  to  rotors  in  compressible  flow 
requires  that  a  physical  model  be  postulated  for  the  planform  of  the  comple¬ 
mentary  wing.  This  planform  is  not  immediately  obvious  as  it  is  for  the 
case  of  wings  in  uniform  flows.  In  applying  the  complementary  wing  concept 
to  calculate  the  effect  of  tip  relief  on  rotors,  full  use  is  made  of  the 
blade  element  theory  while  extending  it,  and  of  the  Prandtl-Glauer t  trans¬ 
formation  with  the  extension  to  rotors  as  has  been  derived  here. 

Blade  element  analysis  of  rotors  is  based  on  the  theory  that  an  elemental 
spanwise  segment  of  a  rotor  blade  acts  as  a  similar  segment  of  a  two-dimen¬ 
sional  airfoil.  The  angle  of  attack  (a) ,  Mach  number  (My)  and  Reynolds 
number  (Re)  of  the  blade  element  are  determined  from  the  motion  at  the  ele¬ 
ment  including  rotation,  advance,  flapping,  etc.,  velocities  and  from  the 
induced  velocity.  Since  ench  element  is  considered  independent  of  the  ad¬ 
joining  elements,  the  force  on  the  element  is  determined  from  the  force 
coefficients  of  the  similar  two-dimensional  airfoil  at  the  angle  of  attack, 
Mach  number,  and  Reynolds  number  calculated.  The  total  thrust  and  torque 
of  the  rotor  are  obtained  by  integrating  the  contribution  of  each  individual 
element  along  the  blade  span  and  the  azimuth  of  the  rotor  disc. 

Consider  a  straight,  rectangular  rotor  blade  of  radius  R  and  a  constant 
chord  cD,  rotating  at  an  angular  velocity  ft  in  compressible  flow,  as  shown 
in  Figure  13  (a).  A  specific  element  is  located  at  a  distance  r  from  the 
center  of  rotation,  i.e.,  a  distance  y  =  R  -  r  from  the  tip,  and  the  ele¬ 
ment  is  at  a  local  Mach  number  denoted  by  My. 
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The  application  of  the  Prandtl-Glauert  type  transformation  to  the  rotor 
in  compressible  flow  results  in  a  nonuniform  chord  blade  in  the  incompres¬ 
sible  flow,  as  was  illustrated  in  Figure  4.  Now  the  complementary  wing 
has  to  be  applied  and  the  question  becomes  that  of  postulating  the  planform 
of  the  complementary  wing.  Attempts  to  apply  a  nonuniforra  chord  comple¬ 
mentary  wing  led  to  physical  and  mathematical  difficulties  and  thus  was 
discarded.  If  instead,  blade  element  theory  is  applied  directly,  then 
each  element  of  the  blade  can  be  considered  as  part  of  an  infinite  span 
airfoil  at  the  local  velocity,  Mach  number,  and  local  chord  in  uniform 
flow.  The  blade  element  theory  may  be  extended,  so  that  instead  of  an 
infinite  span  airfoil,  the  aerodynamics  of  the  rotating  element  is  char¬ 
acterized  by  an  element  located  nt  the  same  spanwise  distance  from  the 
tip  of  a  finite  wing,  whose  span  is  the  same  as  that  of  the  rotating 
blade.  The  complementary  wing  then  is  of  constant  chord  to  infinity  and 
the  results  of  the  previous  chapter  may  be  applied.  To  summarize,  an 
element  at  r  of  a  blade  rotating  in  compressible  flow  at  a  local  Mach  num¬ 
ber  My  has  the  same  tip  relief  effect  as  an  element  located  at  R  -  r  of 
a  similar  wing  of  span  R,  with  chord  cQ  //l  -  My .  This  is  illustrated 
in  Figures  13  (a)  through  (c) . 

This  may  also  be  approached  by  applying  the  blade  element  theory  first 
and  only  then  the  Prandtl-Glauert  transformation.  Since  both  transforma¬ 
tions  are  linear,  their  order  is  immaterial.  Using  the  extended  blade 
element  theory,  an  element  of  the  rotating  blade  has  the  same  aerodynamic 
characteristics  as  an  element  at  a  distance  y  from  the  tip  of  a  wing  with 
span  R,  with  constant  chord  c0,  travelling  at  a  uniform  velocity  Sir  and 
Mach  number  My.  To  this  wing  is  then  the  Prandtl-Glauert  transformation 
applied,  changing  the  constant  chord  to  a  new  value  c0  /  /  1  -  M^.  The 
final  result  i6  physically  and  mathematically  the  same  from  either  approach, 
although  the  latter  may  be  more  satisfying  physically  at  the  intermediate 
step. 

In  both  approaches  the  assumptions  made  in  deriving  the  Sears12  relation 
between  rotating  and  uniform  flow  are  not  thoroughly  observed.  In  the 
first  approach,  a  cylinder  with  a  varying  chord  is  transformed  to  plane 
flow;  while  in  th.  second, the  transformation  is  applied  to  compressible 
flows.  These  extensions  of  Sears'  theoretical  work  are  those  that  are 
used  in  all  blade  element  analyses.  In  addition,  the  transformation  is 
used  for  the  finite  span  case. 

Some  theoretical  validity  for  the  blade  element  approach  as  extended  to  the 
use  with  finite  span  was  recently  given  by  Caradona  and  Isom.56  They  cal¬ 
culated  the  chordwise  pressure  distribution  for  a  nonlifting  rotor  in  com¬ 
pressible,  inviscid  flow  and  compared  it  to  the  pressure  distribution  on 
a  wing  at  a  similar  6panwise  position  at  the  tip  Mach  number.  The  results 
show  good  agreement,  except  in  areas  affected  by  shocks. 

An  indication  of  the  need  for  the  use  of  finite  rather  than  infinite  aero¬ 
dynamic  characteristics  in  blade  element  theory  may  be  found  in  this  paper 
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Figure  13.  Transformation  of  a  Rotor  Blade  Element. 
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by  Caradona  and  Isom.56  The>  compare  the  chordwise  pressure  distribution 
at  similar  spanwise  stations  on  a  rotor  and  on  a  wing  at  Mach  number  of 
the  rotor  tip.  This  differs  from  the  blade  element  theory  extended  here, 
as  a  wing  at  a  different  Mach  number  is  to  be  used  at  each  station.  For 
the  first  one-half  semichord,  there  is  no  perceptible  difference  in  their 
figure  presenting  this  result.*  Further  inboard,  the  rotor  has  a  lower 
(negative)  pressure  coefficient  (i.e.,  smaller  perturbation  velocities) 
than  the  wing  at  the  tip  Mach  number.  This  is  also  true  at  all  stations 
when  the  Mach  number  is  such  as  to  cause  local  supersonic  speeds.  The 
result  may  be  expected  since  the  effective  Mach  number  on  the  blade  is 
lower  than  on  the  wing  to  which  it  is  compared. 


For  compressible  flow,  we  must  further  apply  the  local  Prandtl-Glauert 
type  transformation  derived  previously.  As  is  usual  in  fixed-wing  type 
flows,  this  correction  is  made  to  the  chord: 

Cc  “  Ci  ey  (105) 

thickness  ratio: 


&  r  ■  t. 
y  c  i 


t1/c1 


aspect  ratio: 


perturbation  velocity: 

Au  «  Au./82 
c  i  y 

The  free-stream  velocity  does  not  change: 


=  Xl/  8y 


U 


OOQ 


U 


“i 


(106) 


(107) 


(108) 


(109) 


Therefore,  the  correction  factor  changes: 

G  »  G  ,/62  =  Au  /U  (110) 

X  xi  y  c  00 


TIP  RELIEF  CALCULATION  IN  BLADE  ELEMENT  ANALYSIS 


The  tip  relief  effect  is  included  quantitatively  in  blade  element  analysis 
by  correcting  the  two-dimensional  airfoil  data  in  accordance  with  the  theory 
presented  here.  The  method  for  making  this  application  in  blade  element 
computer  programs  will  be  described  in  a  step-by-step  manner. 


*For  M  .  »>  0.85  the  difference  is  less  than  0.01  at  this  station, 

tip 


1.  In  blade  element  analysis  computer  programs,  the  airfoil  data, 
especially  the  drag  coefficient,  is  available  in  some  form  as 

a  function  of  angle  of  attack,  Reynolds  number,  and  Mach  number. 

2.  At  the  blade  element  under  consideration,  the  velocities  are 
determined;  thus  angle  of  attack,  Reynolds  number,  and  Mach 
number  can  then  be  determined. 

3.  The  correction  factor  for  the  blade  element,  knowing  the  distance 
of  the  element  from  the  tip,  is  calculated  using  the  method 
described  before.  Here,  one  can  represent  the  wing  as  a  semi¬ 
infinite  or  finite  aspect  ratio  wing  and  use  one  or  three  terms 
of  the  Taylor  series  or  a  numerical  integration  scheme  to  obtain 
the  correction  factor. 

4.  Calculate  by  *  1  -  (f:r/a  )2  for  the  element  and  transform  the 
correction  factor  to  the“compress.lble  case  (see  equation  (110)  ). 

5.  Find  the  effective  Mach  number  (Me£f) >  using  equation  (54). 

6.  Obtain  from  the  airfoil  data  the  drag  coefficient  at  the  a  and  Re 
obtained  in  step  2  and  at  the  Mgf£  obtained  in  step  5. 

7.  Calculate  ACD/CD,  using  equation  (58). 

8.  Determine  CDeff,  combining  the  drag  coefficient  found  ir  >  6 

and  the  drag  correction.  ACq/Cq  found  in  step  7.  Use  the  ..fective 
drag  coefficient  in  calculation  of  the  drag  force  of  the  e  .ement. 

These  steps  re,  of  course,  repeated  for  each  segment  each  time  as  they  are 
calculated  ing  the  iterations. 

A  SIMPLIFIED  BLADE  ELEMENT  ANALYSIS 


To  obtain  an  indication  of  the  magnitude  and  variation  of  the  tip  relief 
effect,  a  greatly  simplified  blade  element  analysis  is  derived  for  the 
case  of  a  hovering  rotor.  These  gross  simplifications  allow  the  use  of  a 
short  computer  program  to  make  the  integration  over  the  rotor  radius.  The 
assumptions  to  be  m  will  be  discussed  in  turn  as  they  are  made. 

The  first  assumption  to  he  made  concerns  the  form  of  the  drag  coefficient 
and  its  variation  with  Mach  number.  It  is  often  (for  example.  Reference  55) 
approximated  by  the  equations 


CD  =  CDi 


0<M<M 


crit 


'D  "  "Di  +  ac  '  Mcrit)  Kcrit 


Cn  —  Cr 


<M<1.  0 


(112) 


or  Infinite  aspect  ratio.  For  an  NACA  U(Jl2  alrtoii  section,  moat 
commonly  ”sed  in  rotor  blades, 


M  =  0.8 

crit 

a  =0.57 
c 


(113) 


As  r_ eviously  noted,  the  correction  for  tip  relief  is  made  up  of  two  parts: 
a  change  in  the  Mach  number  AM  and  a  change  in  the  arag  coefficient  AC^, 
as  obtained  from  equations  (54)  and  (58)  respectively.  An  inspection  of 
Figure  9  will  indicate  that  the  change  in  Mach  number  is  by  far  more 
important  than  the  change  in  drag  coefficient.  Note  that  the  change  in 
Mach  number  assures  a  change  in  the  drag  coefficient  in  the  subsonic  drag 
rise  range  of  Mach  numbers.  Therefore,  in  this  simplified  analysis  the 
change  in  drag  coefficient  will  be  neglected. 


From  equation  (112),  the  drag  coefficient  for  a  finite  wing  may  be  written  as 


dfw 

or  rearranging  slightly. 


=  CD  +  a  ((M  -  AM)  -  M  .  )  +  AC 

1  c  crit  1) 


CDFw  =  CDi  +  ac  (  (M  -  Mcrit>  -  M>  +  CI 


(114a) 


(114b) 


Inspection  of  this  equation  indicates  that  the  drag  coefficient  of  a  finite 
wing  is  made  up  of  three  parts: 

1.  The  incompressible  drag  of  an  infinite  wing 

2.  A  correction  for  the  compressible  drag  rise  (M  -  M  ) 

c  crit 

3.  A  correction  for  the  tip  relief  effect  made  up  of  two  terras 

'v  -a  AM  +  ACn 
c  u 

The  ACD  term  may  he  neglected  for  the  reasons  discussed  above.  It  should 
be  pointed  out  that  by  neglecting  the  aC^  terra,  the  tip  relief  effect  is 
assumed  negligible  below  M  ,t,  since  the  drag  coefficient  is  independent 
of  Mach  number  as  stated  in  the  first  of  equation  (112). 


The  torque  of  an  element  of  a  rotor  blade  due  to  drag  is 

dQ  =  rclD  co94> 

-  p  (fir)2  2c  CD  r  dr 


(115) 


The  usaal  small-angle  assumption  is  made.  Since  the  equation  (equations 

(114) ). !or  the  drag  coefficient  is  linear,  upon  substituting  into  equation 

(115) ,  Che  change  in  torque  on  the  element  due  to  tip  relief  may  be 
separated  and  is 
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(116) 


dAQ  =  p(Gr)2  2c  ac  AM  rdr/2 

or  in  coefficient  form  for  b  hovering  rotor  blades 

dAC  =  7  ““  a  AM  r3dr 
Q  2  nR  c 


(117) 


From  the  definition  of  solidity  and  the  transformation  of  equations  (105) 
and  (107) , 


2bcr 

2b 

bB 

rR 

=  n  j  6  =  0  , 

tiE  i  y  y  1 

UA 

(118) 


Using  equations  (54)  and  (118)  in  (117), 


Integrating , 


dACg  7s-  ac(l  +  -^M2)  (Gx/62)  r3  dr 


M2)Gxr3dr 


(119) 


(120) 


It  should  be  emphasized  that  M,  6  and  Gx  are  all  functions  of  the  radiai 
position  along  the  blade.  Also  note  that  Gx  depends  on  F,  the  airfoil 
coordinates,  whose  thickness  must  be  transformed  in  accordance  with 
equation  (106). 


In  equation  (120)  the  correction  factor  needs  to  be  calculated  separately. 
In  the  simplest  case,  the  correction  factor  is  obtained  by  assuming  a  semi¬ 
infinite  wing  and  using  the  first  term  of  the  Taylor  series  expansion.  The 
change  in  Mach  number  then  becomes,  using  equation  (23b)  v35c)  and  (54), 

Y-l  1 


AM/M  =  (1  + 


M2)  2(1- 


:) 


/(c/y)2  +  1 


I1/(4ng2) 


=  2(1  + 


Y-l 


M2) 


/(l/(6  >c(l-f)2)2  +l] 


h_ 

T  4ttS2 


(121) 


The  factor  2  occurs  because  2X=R.  The  value  of  I^/i  is  given  in  Table  II 
Substituting  into  equation  (116)  results  in  the  integral 
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dr 


(122) 


Mr 


(  1  +  M2) 


Even  in  the  simplest  case  of  using  the  correction  factor  as  for  the  first 
term  of  the  Taylor  series  for  a  semi-infinite  wing,  as  shown  in  equation 
(122)  ,  the  integration  of  equation  (120)  appears  to  be  too  complicated  to 
be  analytically  integrated.  Therefore,  three  computer  programs  to 
numerically  carry  out  the  integration  were  written;  they  are  listed  in 
order  of  complexity: 

1  A  program  using  only  the  first  term  of  the  Taylor  series  for  a 
semi-infinite  wing  only  as  the  correction  factor,  as  in  equation 
(122). 

2.  A  program  using  the  Taylor  series  as  the  correction  factor  with 
options  for  using  only  the  first  or  all  three  terms  of  the  series 
and  a  finite  or  a  semi-finite  wirg. 

3.  A  program  using  the  numerically  integrated  version  of  the  correction 
factor. 


These  programs  are  similar  and,  therefore,  will  be  described  together.  They 
use  the  trapezoidal  rule  for  the  integration,  with  a  maximum  of  50  steps 
and  the  values  of  the  Independent  variable  (rotor  radius)  given  as  input 
data.  Each  program  consists  of  a  main  program,  which  sets  up  the 
necessary  variables  and  contains  the  integration  loop.  To  calculate  the 
integrand,  a  function  subroutine  is  called.  The  correction  factor  is 
calculated  within  the  function  in  program  1;  for  the  other  two  programs, 
the  correction  factor  is  calculated  in  a  further  subroutine. 

In  the  third  program,  the  correction  factor  is  calculated  by  numerically 
integrating  equation  (16)  as  described  on  page  15.  For  parabolic  arc 
airfoil  shapes,  22  evenly  divided  segments  were  used  along  the  chord  from 
0.0000  to  0.9999  (since  this  airfoil  is  symmetric  about  the  £  =  0  axis), 
while  30  points  from  +0.3999  to  -0.9999  were  used  for  the  NACA  0012  airfoil. 

Because  of  the  assumptions  made  concerning  the  drag  coefficient,  there  is 
no  effect  on  the  torque  due  to  tip  relief  on  any  segment  where  the  local 
Mach  number  is  below  eight-tenths.  Also,  there  is  a  strong  singularity 
at  tip;  therefore,  the  upper  limit  is  0.999950  rather  than  exactly  one. 
Forty-nine  unevenly  spaced  segments  are  used.  Using  only  39  segments  in 
some  test  cases  resulted  in  a  change  only  in  the  third  significant  figure, 
but  since  there  is  only  a  small  increase  in  computing  time,  49  segments 
are  used.  The  intervals  ware  chosen  so  that  the  values  of  the  integrals 
from  the  segments  are  roughly  equal. 
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RESULTS 


The  crjputer  program  discussed  previously  were  run  with  a  number  of 
sets  of  data,  with  the  various  options  available  for  calculating  the 
correction  factor.  The  results  from  these  computer  runs  are  discussed, 
as  is  the  problem  of  comparison  with  experiment. 

In  making  these  calculations  many  major  assumptions  were  made,  the  most 
important  being  the  simplified  form  of  the  drag  coefficient  and  neglecting 
the  ACd  term  in  the  tip  relief  correction  on  page  56.  These  assumptions 
are  such  that  there  is  no  variation  of  torque  due  to  tip  relief  with  lift. 

To  make  full  use  cf  the  theory  presented  here,  it  should  be  applied  to  the 
drag  coefficient  used  in  a  strip  analysis  computer  program. 

It  should  be  noted  in  regard  to  the  accuracy  required  that  the  calculation 
is  of  the  change  in  torque  coefficient  due  to  tip  relief.  The  order  of 
magnitude  of  Cn  is  about  25  x  10-s ,  while  at  the  same  time  the  change  in 
torque  coefficient  due  to  tip  relief  is  about  2  x  10~s.  Thus,  there  is 
a  factor  of  10  between  these  two  numbers,  which  are  subtracted  from  each 
other.  Three-place  accuracy  on  ACq  requires  four-place  accuracy  of  Cq. 

Results  from  the  computer  output  indicate  the  differences  between  the 
various  means  of  calculating  the  tip  relief.  This  is  summarized  in  Table  IV, 
The  choices  are  the  use  of  finite  and  semi-infinite  wing,  and  the  use  of 
numerical  integration  and  one-  and  three-term  Taylor  series  approximation. 
Further,  calculations  were  made  using  parabolic  arc  and  NACA  four-digit 
series  airfoil. 
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TABLE  IV.  CHANGE  IN  TORQUE  DUE  TO  TIP  RELIEF 

Tip  Mach 

Solidity  Per 

-AC_  Per  Blade 

Q 

Number,  M 

Blade 

Finite  Wing 

Semi-Infinite 

Wing 

Using  three  terms  of  t 

he  Taylor  series 

O.85 

.0159 

1. 026x1 0~5 

1.02l+xl0"5 

0.90 

.0159 

1.8h7 

1.81+1 

0.95 

.0159 

3.812 

3.796 

0.85 

.0232 

1.836 

1.830 

0.90 

.0232 

3.369 

3.352 

0.95 

.0232 

6.7h8 

6.699 

Using  numerical 

integration 

0.85 

.0159 

.993 

.991 

0.90 

.0159 

1.789 

1.781+ 

0.95 

.0159 

3.708 

3.693 

Using  the  first  term  of 

the  Taylor  series 

0.85 

.0232 

1.718 

1.711 

0.90 

.0232 

3.180 

3.162 

0.95 

.0232 

6.203 

6.153 

Parabolic  arc 

airfoil 

The  differences  between  the  use  of  finite  and  semi-infinite  wings  for 
calculating  ACq  are  very  small.  The  differences  are  in  the  third  signifi¬ 
cant  figure,  and  thus  it  is  recommended  that  calculations  be  made  using  the 
semi-infinite  wing  concept.  Further  discussion  will  be  restricted  to  the 
semi-infinite  case. 


Figure  l1*  shows  the  change  in  torque  coefficient  due  to  tip  relief  obtained 
from  the  three  means  of  calculating  the  correction  factor.  The  numerical 
integration  requires  much  more  computing  time,  but  since  the  differences  ere 
fairly  small,  it  is  believed  that  the  use  of  the  first  term  of  the  Taylor 
series  will  give  satisfactory  results. 
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The  change  in  torque  varies  drastically  with  tip  Mach  number,  and  increases 
rapidly  as  Mach  one  is  approached.  This  is  due  to  the  use  of  the  compressi¬ 
ble  flow  transformations  discussed  on  page  50.  All  of  these  transformations 
behave  in  a  manner  similar  to  the  Frandtl-Glauer t  transformation,  having 
an  infinite  singularity  at  Mach  number  equal  to  one.  The  results  presented 
here  used  the  Prandtl-Glauert  transformation,  although  for  application, 
the  extensions  cf  the  Prandtl-ulauer t  transformation  that  were  discussed 
should  be  investigated. 

The  change  in  torque  coefficient  divided  by  solidity  per  blade  due  to  tip 
relief  at  two  different  solidities  is  shown  in  Figure  15.  Solidity  of 
0.0159  per  blade  corresponds  to  a  blade  of  aspect  ratio  20,  while  the 
solidity  per  blade  of  0.0234  represents  the  two-bladed  UH-1F  helicopter. 

The  quantity  ACq/(c/b)  increases  with  solidity.  This  is  as  expected;  with 
the  greater  amount  of  blade  area,  the  flow  is  less  like  two  dimensional. 

Also  indicated  on  Figure  15  is  some  comparison  with  test  data.  The  data 
is  obtained  from  Reference  31,  Figure  9,  where  the  measured  power  coef¬ 
ficient  (numerically  equal  to  the  torque  coefficient)  of  a  UH-1F  heli¬ 
copter  at  u  5  0.25  is  presented  along  with  the  predicted  performance 
based  on  Tanner.23  The  experimental  value  of  the  torque  coefficient  due 
to  tip  relief  is  the  difference  between  these  two  curves. 

The  figure  shows  remarkable  agreement,  but  because  of  the  difference  in 
conditions  used  for  the  two  cases,  care  must  be  taken  in  interpreting  this 
figure.  First,  the  difference  in  scale  between  the  two  figures  is  very 
large,  the  data  being  taken  from  a  small  figure.  Thus,  the  curve  indi¬ 
cating  the  test  should  show  a  wide  error  band.  At  the  same  maximum  tip 
Mach  number,  a  larger  portion  of  a  hovering  rotor  is  in  compressible  flow 

when  compared  to  an  advancing  rotor.  Thus,  the  hovering  rotor  may  be  ex¬ 

pected  to  have  a  larger  tip  relief  correction  than  an  advancing  rotor  at 
the  same  tip  Mach  number,  as  is  indicated  <n  Figure  15.  This  figure  does 
show  that  the  theory  compares  well  with  flight  test  data  in  both  order  of 

magnitude  and  trend  with  Mach  number.  Thus  the  theory  presented  here  can 

help  to  explain  differences  between  prediction  and  test  data  of  helicopter 
rotors  operating  in  the  compressible  flow  regime. 
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CONCLUSIONS 


A  derivation  has  been  presented  which  attempts  to  account  for  differences 
between  performance  prediction  and  flight  teat  data,  as  being  due  to  tip 
relief  effect.  The  theory  is  based  upon  the  complementary  wing  concept 
and  results  in  a  correction  factor  to  the  local  free-stream  velocity. 

Small  perturbation  theory  for  rotors  in  compressible  flow  is  discussed  to 
show  how  a  Prandtl-Glauert  type  transformation  may  be  used  in  rotary-wing 
theory.  Application  of  the  correction  factor  with  the  Prandtl-Glauert  type 
transformation  in  blade  element  analysis  is  discussed.  The  use  of  the 
theory  in  a  simplified  blade  element  analysis  results  in  a  change  in 
torque  due  to  tip  relief,  which  is  compared  with  data. 

Results  from  the  simplified  analysis  are  of  the  same  order  of  magnitude 
as  the  discrepancy  in  test  data,  and  indicate  the  same  trend  with  Increasing 
tip  Mach  number. 

To  fully  exploit  the  theory,  it  should  be  made  a  part  of  a  performance 
computer  program  using  blade  element  analysis.  By  doing  so,  more  meaningful 
comparison  with  flight  test  data  may  be  made. 

The  tip  effect  is  only  one  of  the  factors  which  influence  the  compressible 
flow  over  helicopter  rotors.  Therefore,  the  quantitative  analysis  presented 
in  this  report  cannot  by  itself  lead  to  a  complete  understanding  of  the 
performance  of  helicopter  rotors  in  compressible  flow. 
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